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ABSTRACT 


Experimental  and  theoretical  researches  were  made  on  the  problem 
of  predicting  the  fatigue  life  of  a  resonant  structure  exposed  to  Jet 
noise  by  testing  It  with  Intense  sound  from  a  siren.  One  panel  design 
of  Alclad  2024  aluminum  was  tested  to  fatigue  with  constamt-amplltude 
and  variable  amplitude  siren  sounds  and  with  jet  noise,  A  close  cor¬ 
relation  of  all  three  results  was  found. 

A  comparison  was  also  made  of  fatigue  lifetimes  of  a  resonant 
cantilever  beam,  of  plain  2024  aluminum,  measured  with  constant  exci¬ 
tation  amplitude  (pure  tone)  and  with  random  excitation.  Random  life¬ 
times  were  found  to  be  shorter  than  predictions  from  constant-amplitude 
data. 
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SECTION  I 


INTRODUCTION 


This  report  describes  results  of  experimental  research,  with 
associated  theoretical  investigations,  on  the  problem  of  predicting 
the  fatigue  life  to  be  expected  of  a  resonant  structure  exposed  to 
Jet  noise  by  means  of  the  results  of  tests  with  intense  sound  from 
a  siren. 

Several  different  testing  methods  and  prediction  schemes  were 
tried.  The  actual  predictions  are  compared  in  section  VII  with  life¬ 
times  measured  in  a  Jet  noise  field. 

The  theoretical  bases  of  the  various  prediction  schemes  are  des¬ 
cribed  in  section  II.  One  scheme  depends  upon  standard  fatigue  curves, 
constant-amplitude  fatigue  tests  with  a  siren,  response  measurements, 
and  theoretical  calculations.  A  second  scheme  is  thoroughly  experi¬ 
mental,  being  based  on  variable-amplitude  fatigue  tests  of  the  struc¬ 
ture  with  a  siren  and  various  response  measurements. 

Some  theoretical  analyses  of  the  effect  upon  fatigue  life  of  a 
nonlinearity  of  the  hard-spring  type  are  also  presented  in  section  II. 

Section  III  contains  experimental  comparisons  of  life  under 
constant-aimplltude  and  random  excitations  for  a  resonant  cantilever 
beam  of  2024  aluminum  alloy.  The  results  with  random  excitation  are 
compared  with  predictions  from  the  constant-amplitude  data.  These 
experiments  represent  a  test  in  miniature  of  one  prediction  scheme, 
under  very  carefully  controlled  conditions. 

Results  of  laboratory  studies  of  the  response  of  the  structural 
panels  used  in  the  full-scale  tests  are  given  in  section  IV,  along 
with  some  response  data  from  the  field  tests. 

The  details  of  the  design  for  the  program  used  in  the  varlable- 
sunplltude  fatigue  tests  are  presented  in  section  V. 

The  results  of  fatigue  tests  with  constant-amplitude  and 
variable-amplitude  sound  from  the  siren  and  of  field  tests  with  Jet 
noise  excitation  are  given  in  section  VI. 


Manuscript  released  by  authors  October  196I  for  publication  as  an 
ASD  Technical  Report. 
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The  results  of  this  research  have  Indicated  a  number  of  areas 
which  require  detailed  Investigation  In  order  to  Improve  the  reli¬ 
ability  of  siren  testing  methods.  These  areas  are  described  In 
section  VIII. 


Throughout  this  report  the  amplitude  of  strain  e.  Is  reported  ln_^ 
logarithmic  measure.  This  "strain  level"  Is  defined  as  20  log-Q(e/lO~ 
the  unit  Is  the  "decibel",  abbreviated  to  db,  and  the  reference  level 
Is  1  mlcrostraln,  l.e.,  10"^.  This  procedure  Is  natural  to  persons 
familiar  with  acoustical  measurements  and  was  found  almost  essential. 

In  view  of  the  available  Instrumentation,  In  order  to  avoid  errors. 


) 
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SECTION  II 


SIREN-JET  EQUIVALENCES  FOR  FATIGUE 


In  this  section  we  consider  the  bases  of  various  methods  for 
predicting  the  lifetime  of  a  structure  exposed  to  Jet  noise  from 
results  of  tests  performed  with  the  more-or-less  pure  tone  sound 
from  a  siren.  The  purpose  of  any  of  these  methods  Is  the  prediction 
of  sound  excitations  In  the  two  environments  which  are  "equivalent" 

In  the  sense  that  the  resulting  fatigue  lifetimes  are  equal.  Alter¬ 
natively,  a  satisfactory  "equivalence"  will  be  achieved  If  the  two 
lifetimes  differ  by  a  predictable  eunount. 

Structures  susceptible  to  sonic  fatigue  damage  are  highly  reso¬ 
nant,  characterized  by  one  or  more  resonance  frequencies.  The  forces 
exerted  by  Jet  sound  are  random  In  nature.  Thus  any  equivalence 
scheme  must  take  cognizance  of  the  differences  In  response  to  random 
and  to  pure  tone  forces,  and  also  of  the  differences  In  fatigue  damage 
done  by  these  responses.  Differing  degrees  of  Idealization  of  those 
two  aspects  of  the  equivalence  question  —  response  and  fatigue  — 
account  for  the  differences  between  the  various  schemes. 

There  Is  one  other  Important  criterion  for  a  satisfactory  slren- 
Jet  equivalence  scheme:  It  must  be  capable  of  experimental  realiza¬ 
tion  with  a  reasonable  Investment  of  equipment  and  effort. 

In  this  section  we  describe  a  number  of  different  equivalence 
schemes.  Some  are  based  more  or  less  directly  on  methods  described 
In  the  literature,  while  others  were  developed  on  this  project. 


A.  LINEAR,  ONE-MODE,  NON-INTERACTING  SYSTEM 

This  equivalence  scheme  Is  valid  for  a  linear  system  resonant  In 
a  single  mode  when  there  are  no  fatigue  Interactions  between  peak 
strains  of  different  levels  (l.e..  Miner's  law  of  fatigue  accumula¬ 
tion  Is  valid).*  The  fundamental  analysis  was  carried  out  by  Miles 
( ref .  1 ) . 

1.  Analysis  for  Fatigue  Daunage 

We  assume  that  the  Jet  sound  pressures  are  random  and  normally 
distributed  with  a  smooth  spectral  density  In  the  narrow  bandwidth  of 
the  structure's  resonance.  These  have  been  demonstrated  to  be  good 
assumptions  (ref.  2).  Then  the  probability  density  of  the  peaks  of 
response  strain  e  at  any  point  Is  a  Rayleigh  distribution  curve 


♦Throughout  this  report  we  use  strain  rather  than  stress  as  the 
measure  of  response  for  reasons  that  are  discussed  below. 
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p(x)  =  X  exp  {-x^/2)  (1) 

where  x  Is  the  ratio  of  the  magnitude  of  a  strain  peak,  e,  to  the  rms 
strain  a: 

X  =  €/a  .  (2) 

The  strain  cycles  are  fully  reversed,  except  when  the  peak  strain  Is 
very  small,  and  the  average  Interval  between  peaks  Is  the  period  of 
natural  vibration. 

In  addition  to  the  assumption  of  no  fatigue  Interactions,  let  us 
assume  that  the  constant-amplitude  fatigue  curve  (peak  strain  versus 
life)  Is  a  straight  line  In  a  log-log  plot: 

=  constant,  (3) 

where  N  Is  the  fatigue  lifetime  In  cycles.  Then,  as  Miles  has  shown, 
particularly  simple  results  are  found.  The  average  damage  per  cycle 
of  random  response  Is  the  same  as  the  damage  per  cycle  of  an  equiva¬ 
lent  response  whose  peak  value  Is 

peak  =  Vs  ((a/2)!l^/“  ,  {4a) 

and  whose  rms  value  Is 

rms  =  ((a/2)!]^/“  .  (4b) 

Steady  sinusoidal  excitation  to  this  equivalent  response  level,  with 
fully-reversed  cycles,  should  therefore  yield  the  same  lifetime  as 
the  random  excitation. 

In  view  of  the  curvature  that  appears  In  experimental  constant- 
amplitude  fatigue  curves.  It  Is  reasonable  to  modify  this  procedure 
by  substituting  the  experimental  curve  for  Eq.  (3).  By  Miner's  law, 
the  fractional  damage  per  cycle  of  peak  strain  e  is  [l/M(e)].  There¬ 
fore  the  average  damage  per  cycle  of  random  response  Is  the  Integral 

00 

1/Nr  =  /  [1/K(e) ]  p(x)  dx  (5) 

o 

and  equals  the  reciprocal  of  the  expected  random  lifetime  In  cycles, 

Nj,.  This  procedure  has  been  developed,  discussed,  and  applied  by  a 
number  of  researchers  (refs.  3>  5)»  The  result  of  the  Integration 

la  a  value  of  random  lifetime  Ng  for  each  value  of  rms  response  strain  a. 
Specific  examples  of  such  calculations  will  appear  In  later  sections  of 
this  report-  (Figs.  8  and  21  below.) 
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2.  Stress-Life  or  Strain-Life? 


In  actuality,  previous  researchers  have  uniformly  used  a  stress- 
lifetime  fatigue  curve  In  these  computations,  whereas  we  recommend 
the  use  of  strain-lifetime  curves.  There  Is,  of  course,  no  difference 
as  long  as  the  strain  In  the  fatiguing  area  Is  below  the  yield  point. 
However,  beyond  the  yield  point,  the  peak  stress  and  peak  strain  are 
nonllnearly  related;  an  essential  conflict  arises  with  the  assumption 
of  linearity  of  response,  made  In  the  preceding  analysis.  The  ques¬ 
tion  to  be  answered  Is:  Which,  If  either,  of  the  two  quantities, 
local  stress  and  local  strain,  can  be  considered  linearly  related  to 
the  overall  structural  response? 

The  answer  depends  on  the  structure  Involved.  In  most  cases  of 
sonic  fatigue,  the  area  of  large  strains  Is  very  localized.  First, 
the  structure  Itself  Is  a  sort  of  strain-concentrating  mechanism  In 
the  sense  that  the  gross  strain,  such  as  might  be  measured  by  a  strain 
gauge,  varies  from  area  to  area  on  It.  Secondly,  the  fatigue  damage 
usually  occurs  In  a  very  localized  region  of  further  strain  concentra¬ 
tion,  at  the  edge  of  a  hole  or  In  a  sharp  bend.  In  these  circumstances. 
It  appears  reasonable  to  assume,  as  a  first  approximation,  that  yielding 
In  the  small  fatiguing  area  does  not  affect  the  surrounding  material  and 
that  the  local  strains  In  the  fatiguing  area  are  linearly  related  to  the 
response.*  We  note  that  Gunn  has  successfully  used  a  similar  approxima¬ 
tion  In  his  analysis  of  fatigue  of  specimens  with  a  stress  concentration 
(ref.  6). 

The  distinction  between  stress  and  strain  becomes  of  particular  Im¬ 
portance  at  short  lifetimes  when  a  significant  proportion  of  damage  Is 
due  to  strains  above  the  local  yield  point.  Whereas  a  constant-amplitude 
stress-lifetime  curve  rapidly  becomes  flat  In  this  region,  the  corres¬ 
ponding  strain-lifetime  curve  continues  to  rise  smoothly  toward  the 
ultimate  strain  (the  strain  for  failure  In  a  "static"  tensile  test) 

(ref.  7).  Because  of  this  difference,  computations  of  a  random  fatigue 
curve  based  on  Eq.  (5)  and  constant-amplitude  stress-life  curves  should 
yield  a  predicted  life  that  Is  very  much  too  small  at  the  short  lifetime 
end.  Just  such  a  difference  between  predictions  and  measurements  has 
been  found  In  some  experiments  (ref.  5)* 


*In  the  language  appropriate  to  Impedance  concepts,  the  structure  as  a 
whole  may  be  considered  a  "source",  delivering  strains  to  the  "load", 
l.e.,  the  local  fatiguing  area.  We  believe  that  the  source  Impedance 
la  large,  the  source  acts  as  a  constant-strain  source,  and  that  the 
strains  at  the  Interface  are  not  significantly  affected  by  the  decrease 
In  load  Impedance  which  local  yielding  represents. 
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3.  Fracture  by  Rare,  Large  Response  Peaks 


The  distinction  between  stress  and  strain  Is  also  Important  In 
settling  the  questions  of  damage  due  to  rare,  very  large  responses. 

In  random  response,  response  levels  of  Indefinitely  large  values  must 
be  expected,  although  only  with  very  small  probability.  A  response 
exceeding  the  ultimate  strength  of  the  material  will  cause  Immediate 
fracture.  How  does  one  assess  the  probability  of  such  damage  by 
fracture,  as  contrasted  with  cumulative  fatigue  damage?  The  procedure 
we  Illustrate  here  Is,  admittedly,  rather  rough,  but  should  give  ap¬ 
proximately  correct  answers.  Any  more  exact  assessment.  Incorporating 
all  the  nonllnearltles  of  material  and  response,  would  be  extremely 
complicated. 

Consider  a  resonant  structure  of  2024  alinnlnum  falling  without 
the  benefit  of  stress  concentrations.  Computations  can  be  made  by 
means  of  Eq.  (5)  using  the  strain-life  failure  curves  of  Freudenthal 
(ref.  8),  and  D' Amato  (ref.  7).  For  an  rms  strain  of  1.9x10"^  (rms 
stress  of  20,000  psl)  one  computes  an  expected  lifetime  of  1.6xl05 
cycles  If  he  Ignores  the  rare,  very  high  strains.  However,  In  this 
case  the  ultimate  strength  (64,000  psl)  Is  only  3.2  times  the  rms 
stress.  If  the  Rayleigh  distribution  [Eq.  (l)]  were  valid  for  stress, 
cycles  of  peak  stress  exceeding  the  ultimate  strength  would  be  expected 
to  occur  6  times  In  every  1000  cycles.*  Clearly,  Immediate  fracture 
from  such  "rare"  response  levels  would  be  much  more  likely  than  failure 
by  the  gradual  accumulation  of  fatigue  dsimage.**  The  lifetime  of  the 
structure  would  not  be  determined  by  the  rate  of  accumulation  of  fa¬ 
tigue  damage  but  rather  by  the  rate  of  occurrence  of  large  responses. 

In  the  example  cited  above,  measured  lifetimes  for  an  ensemble  of 
Identical  structures  would  show  a  broad,  flat  distribution  (equal 
probability  density  in  lifetime)  from  zero  to  about  3OO  cycles  (half 
the  average  Interval  between  peaks  exceeding  the  ultimate  strength), 
with  larger  lifetime  becoming  gradually  less  probable. 


p 

*The  probability  of  exceedances  of  a  response  value  x  is  exp[-x  /2]. 

**If  cycles  of  large  response  levels  occur  In  clumps,  the  probability 
of  occurrence  of  a  clump  Is  less  than  the  probability  of  exceedances 
by  a  factor  which  Is,  as  a  first  approximation,  the  average  number  of 
cycles  In  a  clump.  The  average  clump  size  has  been  determined  by  anal¬ 
ysis,  for  both  linear  and  nonlinear  systems  (ref.  9).  For  practical 
values  of  damping,  there  Is  no  significant  tendency  to  clumping  of  peaks 
exceeding  3  times  the  rms  response. 
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The  fact,  demonstrated  In  the  literature  (r^fs.  3>  5)  and  In 
this  project,  that  random  fatigue  tests  of  resonant  structures  can 
be  made  at  short  lifetimes  without  any  suggestions  In  the  results 
of  failure  by  rare  exceedances  of  the  ultimate  stress.  Is  positive 
evidence  of  the  falsity  of  using  stress  as  the  Independent  variable 
In  Eqs.  (l)  and  (5) . 

If  the  Independent  variable  In  Eqs.  (l)  and  (5)  Is  taken  to  be 
strain,  the  problem  of  failure  by  rare  exceedances  of  the  ultimate 
strength  does  not  arise.  In  the  numerical  example  used  above,  the 
strain  resulting  In  fracture  on  the  first  cycle  {about  O.I9  for  2024 
alvunlnum)  so  far  exceeds  the  rms  strain,  1.9xlO~3^  that  a  computed 
probability  of  exceedance  (l  In  lO^^^l)  is  undoubtedly  meaningless. 

In  this  view,  the  only  significant  failure  mechanism  Is  accumulative 
fatigue. 

4.  Analysis  for  Response 

The  preceding  discussion  (subsection  l)  has  outlined  the  evalua¬ 
tion  of  a  lifetime  for  random  fatigue  from  (l)  the  rms  value  of  re¬ 
sponse  strain  and  (2)  failure  curves  (peadc  strain-lifetime)  for  fully 
reversed  strain  cycles  of  constant  amplitude.  A  complete  equivalence 
scheme  must  also  Include  a  means  of  relating  the  rms  response  to  the 
sound  pressures  exciting  the  structure. 

We  shall  not  consider  methods  for  solving  the  response  problem 
by  pure  computation.  To  some  extent  this  Is  due  to  the  natural  bias 
of  an  experimenter;  such  a  procedure  seems  only  to  add  one  more  hazard 
to  a  long  chain  of  hazardous  predictions.  In  any  case,  a  number  of 
reasonable  experimental  methods  for  solving  the  response  problem  are 
available.  (Computational  methods  have  been  adequately  treated  else¬ 
where  ( ref .  10). ) 

In  the  Immediately  succeeding  paragraphs  we  summarize  the  equa¬ 
tions  governing  the  one-mode  linear  response  of  a  structure  to  sonic 
excitation.  Both  pure  tone  and  random  excitations  are  considered. 

The  dynamical  equation  for  one-mode  response  of  a  structure  to 
pure  tone  sonic  excitation  can  be  written  (ref.  11 ) 

^^mech  “  ^tot  “  ^tot 


where  v  Is  the  response  velocity. 


Z 


mech 


Is 


the  vacuo  mechanical  Impedance, 


Ftot  generalized  force, 

T  Is  a  transfer  function,  with  units  of  area, 
Ptot  sound  pressure  acting  on  the  face 


and 

of  the  structure. 
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If  the  structure  is  large  or  the  motion  complicated^  considera¬ 
tion  must  be  given  to  the  dependence  of  T  upon  the  spatial  character 
Istlcs  of  the  sound  field  (ref.  12).  We  shall  assume  here  that  the 
structure  is  small  and  that  the  sound  pressure  is  sensibly  constant 
over  the  face  of  the  structure,  so  that  T  is  therefore  a  constant  in 
all  environments. 

The  total  sound  pressure  can,  by  superposition,  be  decomposed 
into  the  sum  of  the  sound  pressure  measured  when  the  structure  is 
blocked  (kept  from  vibrating)  and  the  radiation  sound  pressure, 
measured  with  the  structure  vibrating  but  with  the  external  sound 
source  turned  off: 


Ptot  =  Pbi  +  Prad^ 


^tot  =  +  ^rad 


(7) 


The  radiation  pressure  is  related  to  the  radiation  Impedance: 

^rad  “  ^rad  “  ”  '^^rad  ^ 
so  that  Eq.  (6)  can  also  be  written 

'^^^mech'^^rad^  ^  '^^tot  ^  ^b£  ^ 

The  sum  of  mechanical  (i^  vacuo)  Impedance  and  radiation  impedance  is 
the  total  impedance  for  vibration  in  the  acoustical  environment.  The 
radiation  Impedance  may  change  when  the  environment  is  changed. 

The  two  dynamical  equations,  Eqs.  (6)  and  (9)^  are  perfectly 
equivalent,  differing  only  in  the  treatment  of  the  radiation  pressure. 
In  Eq.  (6)  the  radiation  pressure  appears  on  the  right  side;  in  Eq.  (9) 
it  appears  on  the  left  in  the  guise  of  radiation  impedance.  The  two 
equations  have  the  same  general  form: 

vZ  =  Tp  .  (lOa) 

For  the  one-mode  case  we  can  write  the  Impedance  as 

Z  =  R  +  j(a)M  -  K/o))  ,  (lOb) 

where  R,  M,  and  K  stand  for  modal  resistance,  mass,  and  stiffness,  and 
where  o)  is  the  frequency  (radians  per  second).  The  response  strain  e 
is  proportional  to  displacement,  so  that 

e  =  C  v/o)  ,  (11) 

where  C  is  an  appropriate  constant  depending  on  mode  shape. 
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The  re^onse  to  pure  tone  excitation  at  the  resonance  frequency 
0)^  =  Is  the  well-known  solution  of  Eq.  (lO): 

€  =  c  v/o)^  =  CTp/o^R,  or 

=  p^CCT/oDqR)^  (12) 

where  we  uSe  the  superscript  bar  to  represent  the  time  average. 

The  mean  square  response  to  random  excitation  is  an  equally  well- 
known  solution  of  Eq.  (lO)  in  the  special  case  where  the  spectral  den¬ 
sity  of  p  is  essentially  constant  within  the  bandwidth  of  response 
(Ref.  l).  That  solution  is 

=  Sp  C^T^/4RK  (13) 

where  S  is  the  spectral  density  of  the  pressure  p  (units  of  pressure- 
squared^per  cps ) . 

We  now  turn  to  various  schemes  for  predicting  the  rms  response  of 
a  structure  in  a  realistic  Jet  environment  from  measurements  in  an 
artificial  environment. 

a.  Low-Level  Random  Excitation 

One  experimental  procedure  for  predicting  the  response  involves 
the  excitation  of  the  actual  structure  by  a  moderately  low  level  of 
random  noise  and  measuring  the  corresponding  low  level  of  response. 
Then  if,  and  only  if,  the  test  environment  is  essentially  identical 
to  the  realistic  environment,  the  ratio  of  rms  response  to  rms  exci¬ 
tation  will  be  constant  in  a  linear  structure.  If  the  two  environ¬ 
ments  differ,  a  suitable  correction  must  be  made.  We  investigate 
here  the  nature  of  the  correction  for  changes  in  the  acoustic  environ¬ 
ment,  that  is  in  the  radiation  impedance. 

(l)  Measuring  blocked  pressure;  If  the  excitation  is  character¬ 
ized  by  the  blocked  sound  pressure,  the  proper  dynamical  equation  is 
Eq.  (9l.  The  prediction  equation  is  then  [see  Eq.  (13)] 

5^)  -  "pjaV  (14a) 

'"(b)  Sibj'^ajV) 

where  R  and  K  are  total  resistance  and  total  stiffness  and  the  sub¬ 
scripts  (a)  and  (h)  refer  to  the  test  and  the  realistic  environments. 
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In  normal  circumstances,  the  total  stiffness  and  total  mass  are  neg¬ 
ligibly  affected  by  the  radiation  impedance.  (This  would  not  neces¬ 
sarily  be  true  If  the  test  environment  Is  an  acoustical  resonance 
chamber.)  In  such  cases  we  may  approximate  Eq.  (l4a)  by 


(l4b) 


where  t|  Is  the  total  loss  factor,  such  as  can  be  determined  from  the 
rate  of  decay  of  vibrations  In  the  absence  of  external  sound.  (The 
quantity  TrT|  Is  the  logarithmic  decrement  In  nepers  per  cycle.) 

(2)  Measuring  total  pressure;  If  the  excitation  Is  characterized 
by  the  total  sound  pressure  at  the  face  of  the  vibrating  structure,  the 
proper  dynamical  equation  Is  Eq.  (6).  A  different  approach  Is  In  order. 

The  first  problem  that  arises  Is  that  the  total  sound  pressure  does 
not  necessarily  have  a  flat  spectral  density;  a  flat  spectral  density  Is 
required  In  order  to  use  the  simple  expression,  Eq.  (13)^  for  mean  square 
response.  Only  If  the  radiation  pressure  Is  essentially  negligible  Is 
one  Justified  In  using  Eq.  (13).  We  now  consider  one  typical  case  In 
more  detail. 


In  a  common  case  of  an  aircraft  panel  vibrating  In  Its  fundamental 
mode,  the  radiation  Impedance  Includes  a  resistance  which  Is  small  com¬ 
pared  with  the  mechanical  (internal)  resistance,  and  a  reactance  corres¬ 
ponding  to  a  virtual  mass  which  Is  small  compared  with  the  mass  of  the 
structure.  However  the  radiation  reactance  Is,  In  the  common  case, 
large  compared  with  the  radiation  resistance. 


When  the  excitation  Is  random,  the  previous  dynaunlcal  Impedance 
equations  go  over  Into  equations  between  the  spectral  densities  of  the 
variables  (Ref.  13).  Thus,  from  Eq.  (9)  we  have  for  the  spectral  den¬ 
sity  of  blocked  pressure 


S 


P 


(bi) 


■'tot 


where  S  Is  the  spectral  density  of  response  velocity.  At  the  reso¬ 
nance  frequency  co  ,  the  total  Impedance  equals  the  total  resistance 
and 


(bi) 


(%)  =  S^(coJr2  /T^ 


(15a) 
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Correspondingly,  from  Eq.  (8),  we  have  for  the  spectral  density  of 
radiation  pressure 


(15b) 


where  In  writing  the  last  line  we  neglect  the  radiation  resistance 
as  small  compared  with  radiation  reactance. 

Now  let  us  look  at  the  ratio  of  Eqs.  (I5a)  and  (l5b): 


where  we  write  the  loss  factor  q  for  R/cuM.  This  ratio,  Eq.  (l5c)  Is 
not  necessarily  small.  Typical  values  for  loss  factor  are 


T|  »  0.02  to  0.04  , 


while  typical  magnitudes  of  the  mass  ratio  for  thin  panels  have  been 
estimated  to  be  (Ref.  l4) 


”rad/”tot 


«  0.01  to  0.1 


Thus  we  conclude  that  the  spectral  density  of  radiation  pressure  may 
easily  equal  or  exceed  the  spectral  density  of  the  blocked  pressure, 
at  the  resonance  frequency. 

Note  that  the  spectrum  of  the  radiation  pressure  has  essentially 
the  same  shape  as  the  spectrum  of  the  response  velocity.  [See  Eq. 
(15b)  and  note  that  the  radiation  Impedance  usually  varies  only 
slowly  with  frequency. ]  Thus  the  energy  of  radiation  pressure  Is 
essentially  limited  to  a  small  frequency  band,  A=ti.  .  f  cps.  On  the 
other  hand  the  spectral  density  of  sound  pressure  is^seidom  determined 
with  a  filter  having  such  a  narrow  bandwidth.  Consider  the  relative 
contributions  of  radiation  and  blocked  pressures  to  the  total  rms 
pressure  measured  In  a  bandwidth  A'=k  f  : 
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If  measurements  are  made  with  filters  common  In  acoustical  research, 
we  have; 

octave-band  filter:  k  =  0.707 

third-octave  filter;  k  =  0.236 

These  values  are  so  much  larger  than  reasonably  expected  loss  factors 
that  the  ratio  of  mean-square  pressures  In  Eq.  (l6}  Is  usually  a  small 
fraction. 

From  this  analysis  we  conclude  that  the  spectral  density  of  radia¬ 
tion  pressure  at  the  resonance  frequency  Is  not  necessarily  negligible 
compared  with  the  spectral  density  of  blocked  pressure.  Even  If  the 
blocked  pressure  has  a  flat  spectrum,  the  usual  case,  the  total  pres¬ 
sure  may  well  not.  Analyses  such  as  that  of  Miles  (ref.  l)  may  there¬ 
fore  not  be  strictly  pertinent.  On  the  other  hand,  the  "spectral 
density"  Inferred  from  measurements  of  total  sound  pressure  In  a  band 
much  wider  than  the  resonance  band  of  the  structure  Is  often  a  good 
estimate  of  the  true  spectral  density  of  blocked  pressure.  Its  use 
In  analyses  such  as  Miles',  using  equations  such  as  Eqs.  (13)  or  (l4), 
can  thereby  often  be  Justified. 

However,  If  In  either  the  test  or  the  realistic  environments  the 
radiation  Impedance  (particularly  the  reactance)  Is  abnormally  high, 
a  more  detailed  analysis  may  be  required. 

b.  Low-Level  Pure  Tone  Excitation 

In  many  cases,  determination  of  the  rms  response  of  measurements 
using  low-level  excitation  with  a  band  of  random  noise  Is  not  conven¬ 
ient.  Then  low-level  pure  tone  sound,  adjusted  to  the  resonance  fre¬ 
quency  of  the  structure,  may  be  used  In  the  test  environment.  The  rms 
response  In  the  realistic  environment  (jet  noise)  can  then  be  predicted 
from  a  combination  of  Eqs.  (12)  and  (13)5 


%  ^t 


r  r 


(17a) 


where  the  subscripts  r  and  t  are  used  to  distinguish  the  "random",  or 
realistic,  and  "tonal",  or  test,  environments. 

We  now  assume  that  the  resonance  frequencies  and  the  modal  masses 
are  essentially  the  same  In  each  environment.  (This  Is  generally 
found  to  be  a  good  assumption  except  when  the  tonal  test  environment 
Involves  a  resonant  acoustic  system.)  Then  Eq.  (l7a)  can  be  rewritten; 
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(17b) 


In  the  special  case  where  the  loss  factor  t)  is  the  same  in  both  en¬ 
vironments,  this  expression  becomes 


r 


(17c) 


a  form  equivalent  to  results  appearing  widely  in  the  literature 
(Refs.  1,  14,  15). 

(l)  To  measure  blocked  or  total  pressure?  The  same  question 
arises  in  this  case  of  tonal  excitation  as  did  in  the  previous  case 
of  random  excitation:  Should  one  measure  the  blocked  sound  pressure 
or  the  total  pressure  (structure  vibrating)?  The  answer  is  somewhat 
different  in  this  case. 

Prom  Eq.  (12)  governing  response  to  a  tone,  it  is  evident  that 
the  ratio 


is  the  same  in  either  case;  the  same  prediction  of  response  will  be 
made.  However,  if  blocked  pressure  is  measured,  the  sum  of  mechanical 
and  radiation  damping  must  be  used  in  the  calculation.  [See  Eq.  (9).] 

If  total  pressure  is  measured,  then  the  mechanical  (in  vacuo)  damping 
alone  must  be  used  in  the  calculation.  [See  Eq.  (6)77  The  difficulty 
with  the  latter  procedure  is  that  the  in  vacuo  damping  cannot  be  meas¬ 
ured;  the  damping  determined  in  the  teaF  environment  is  the  sum  of 
mechanical  and  radiation  damping.  (Any  of  the  standard  techniques  may 
be  used.) 

An  important  simplification  occurs  if  it  can  be  assumed  that; 

1)  the  radiation  damping  in  the  random  environment  is  negligible; 

2)  the  mechanical  damping  is  the  same  in  both  random  and  tonal 
environments. 

Then  the  total  pressure  technique  may  be  used  with  the  simpler  Eq.  (l7c) 
to  predict  the  random  response,  even  when  the  radiation  damping  is  large 
in  the  tonal  test  environment.  In  the  same  situation,  the  blocked  pres¬ 
sure  technique  requires  that  Eq.  (I7b)  be  used  in  order  properly  to 
account  for  radiation  damping. 
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In  contrast  to  this  simplification  resulting  from  the  use  of  total 
pressure,  there  Is  some  question  about  experimental  precision  arising 
from  the  use  of  this  technique.  In  some  cases  the  radiation  damping 
dominates  Internal  mechanical  damping  In  the  test  environment;  this 
occurred  In  the  siren  experiments  reported  below.  Then  the  total  pres¬ 
sure  Is  much  smaller  than  the  blocked  pressure.  The  total  pressure  Is 
the  vectorial  sum  of  the  blocked  pressure  and  the  radiation  pressure. 

At  resonance  all  three  components  are  In  phase  (or  anti-phase)  as  shown 
In  Fig.  1.  The  blocked  pressure  Is  essentially  constant,  for  small 
changes  In  frequency  near  resonance.  However  the  radiation  pressure, 
being  proportional  to  response  velocity,  varies  both  In  phase  and 
magnitude;  Its  trace  In  the  vector  plot  (Fig.  l)  lies  on  a  circle  pass¬ 
ing  through  the  end  of  the  blocked  pressure  vector. 

Now,  the  response  strain  Is  proportional  to  the  response  velocity 
and,  therefore,  to  the  radiation  pressure.  It  appears  from  Fig.  1  that 
the  magnitude  of  the  ratio  ^vP^,  which  appears  In  the  prediction  for¬ 
mula,  Eq.  (I7)i  may  be  very  sensitive  to  tuning  In  some  circumstances; 
the  fractional  change  In  this  ratio  may  be  much  larger  than  the  frac¬ 
tional  change  In  response.  Therefore  a  slight  detuning,  not  noticeable 
In  the  tuning  procedure,  may  lead  to  larger  errors  In  the  prediction 
when  the  total  pressure  technique  Is  used.  This  effect  should  be  In¬ 
vestigated  qualitatively.  (Note  that  tuning  with  phase  measurements 
would  minimize  this  source  of  error. )  We  believe  this  effect  was 
responsible  for  some  of  the  scatter  observed  In  the  measurements  re¬ 
ported  below. 

To  conclude  this  discussion,  either  the  blocked  pressure  technique 
or  the  total  pressure  technique  can  be  used  In  the  prediction.  If  one 
Is  sufficiently  careful  to  verify  the  assumptions.  The  choice  Is  In¬ 
timately  related  to  the  problems  of  determining  the  damping  factor  and 
of  the  variation  In  damping  with  acoustic  environment. 

5.  Summary  of  Equivalence  Scheme 

The  equivalence  scheme  described  above  Is  based  on  an  analysis 
which  assumes  both  a  linear,  one-mode  structure  and  fatigue  accumula¬ 
tion  without  Interactions  between  the  damage  done  by  successive  peaks 
of  strain.  The  procedure  for  applying  this  scheme  can  be  summarized 
Into  three  steps. 

Step  1.  Pressure-Nominal  Strain:  The  rms  value  of  nominal  strain 
(as  measured  by  a  strain  gauge  at  any  convenient  spot  on  the  structure) 
which  will  result  from  excitation  by  a  specified  random  (jet)  noise 
level  Is  predicted  by  response  measurements.  These  tests  can  use  a 
pure  tone  exciting  the  structure  at  resonance.  In  which  case  Eq.  (l?) 

Is  used  for  the  prediction.  Alternatively,  these  tests  can  use  a  ran¬ 
dom  sound.  In  which  case  Eq.  (l4)  Is  used  for  the  prediction.  Some 
care  Is  required  to  get  an  accurate  prediction  If  the  acoustical  en¬ 
vironment  of  the  structure  In  these  response  tests  Is  not  realistic. 
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AT  RESONANCE 


TOTAL  PRESSURE 


BELOW  RESONANCE 

(jJ  <  (jJq 


ABOVE  RESONANCE 


FIG.I  PHASE  AND  MAGNITUDE  RELATIONS  BETWEEN 

BLOCKED,  RADIATION,  AND  TOTAL  SOUND  PRESSURE 
COMPONENTS  FOR  PURE-TONE  TESTING  NEAR 
RESONANCE 
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step  2.  Nominal  Strain  -  Fatigue  Life  (Constant  Amplitude); 

The  structure  is  excited  by  a  siren  to  a  constant  amplitude  of  alter- 
nating  strain,  and  the  fatigue  life  Is  detemined.  From  this  data,  a 
constant  amplitude  fatigue  curve  is  constructed.  One  may  use  fatigue 
curves  derived  for  the  material  by  measurements  on  small  specimens, 
assuming  that  the  strain  concentration  factor  (between  nominal  strain 
and  fatiguing  strain)  is  a  constant.  Alternatively,  one  may  excite  to 
failure  many  nominally  identical  structures  and  construct  a  fatigue 
curve  directly. 

Step  3.  RMS  Random  Nominal  Strain  -  Fatigue  Life;  The  last  step 
is  the  computation  of  a  "random  S-N"  curve  in  which  the  ordinate  is  the 
rms  value  of  the  nominal  strain.  The  computation  is  based  on  the 
assumed  Rayleigh  distribution  of  strain  peaks  [Eq.  (l)]  and  the  inte¬ 
gral  formulation  of  Miner's  law,  Eq.  (5). 


B.  THEORETICAL  EQUIVALENCE  FOR  NONLINEAR  ONE-MODE,  NON-INTERACTING 
SYSTEM 

An  analytical  investigation  was  carried  out,  with  the  aid  of  a 
digital  computer,  to  determine  the  effect  of  nonlinearities  in  the 
response  of  a  structure  upon  the  fatigue  life  (ref.  I6).  The  presence 
of  nonlinearity  leads  to  a  change  in  the  probability  density  of  strain 
peedcs  in  the  response  to  random  excitation.  These  changes  can  lead  to 
a  rate  of  accumulation  of  damage  which  is  different  than  that  for  the 
linear  structure,  even  though  the  excitation  is  the  same.  We  summarize 
that  Investigation  here. 

The  analysis  maintains  the  assumption  that  only  one  mode  of  re¬ 
sponse  is  excited.  Neither  is  any  attempt  made  to  correct  for  fatigue 
Interactions  between  different  strain  levels,  except  to  the  extent  that 
the  correction  can  be  properly  represented  by  a  "fictitious  S-N  curve", 
as  suggested  by  Freudenthal  (ref.  I7). 

The  analysis  tedces  as  a  model  for  the  nonlinear  structure  a  beam 
with  pinned  ends,  responding  in  its  fundamental  mode.  The  displacement 
s(t)  at  the  center  of  such  a  beam  satisfies  the  nonlinear  equation  of  a 
hard-spring  resonator; 

Ms  +  Rs  +  Ks(l+s^/£^)  =  F(t)  ,  (18) 

where  mass  M,  resistance  R,  low-level  stiffness  K,  and  the  "character¬ 
istic  displacement"  I  are  constants.  The  quantity  F(t)  is  the  general¬ 
ized  force.  As  long  as  s  is  small  compared  with  I,  this  is  the  linear 
one-mode  equation.  The  characteristic  displacement  I  equals  twice  the 
radius  of  gyration  of  the  beam's  cross  section  If  the  ends  of  the  beam 
are  firmly  pinned.  Thus,  in  a  flat  plate  of  thickness  t,  £  =  t/VT* 
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The  surface  strain  £  In  the  center  of  the  beam  (the  point  of  majci- 
mum  strain)  is  nonlinearly  related  to  the  displacement  s  by  the  equation 

e  =  cs(1+s/2V3  !>)  ,  (19) 

where  c  is  a  constant  of  proportionality.  The  linear  term  in  Eq.  (19) 
represents  the  usual  bending  strain  which  results  from  a  transverse  dis¬ 
placement  of  the  beam's  center.  The  nonlinear  term  represents  the 
membrane  strain  due  to  the  elongation  of  the  beeim  when  the  transverse 
displacement  is  large.  ► 

1.  Average  Damage  Per  Cycle 

The  average  fatigue  deimage  per  cycle  was  computed  by  Miner's  law 
with  the  further  assumption  that  the  constant-eimplltude  fatigue  curve 
is  a  straight  line,  in  a  log-log  plot,  with  a  slope  (-a).  These  are 
the  assumptions  used  to  derive  Eqs.  (3)  and  (4)  above.  However  the 
exact  probability  density  of  peadcs  was  used  Instead  of  the  Rayleigh 
distribution,  Eq.  (l),  which  holds  for  the  linear  structure. 

The  results  of  these  computations  are  presented  in  Pig.  2.  The 
ordinate  in  Pig.  2  is  the  ratio  of  the  average  daimage  per  cycle  in  non¬ 
linear  random  response  to  the  "linear  estimate"  of  average  deunajge  per 
cycle.  The  linear  estimate  is  the  daimage  that  would  be  sustained  by 
the  same  structure  if  the  characteristic  displacement  I  were  very  large; 
it  is  the  value  that  would  be  predicted  by  the  procedure  outlined  above 
in  subsection  A  [Eqs.  (4)  and  (13)].  The  abscissa  in  Pig.  2  is  a  lin¬ 
earity  parauneter  y  defined  by 

y  =  2^/0^  (20) 

where  a  is  the  linear  estimate  of  the  rms  of  the  displacement  s. 

Linear  response  is  represented  by  y  =  ». 

The  three  top  curves  in  Pig.  2,  labelled  CASE  B,  are  the  results 
for  the  plnned-beaim,  for  values  of  a  =  4,  8,  and  l6.  The  curves  show 
that  the  damage  rate  Increases  with  increasing  nonlinearity  (smaller  y) 
when  the  nonlinearity  is  small.  There  is  a  maximum  point  in  each  case. 

Thus  for  a  value  a  =  l6  (appropriate  to  2024  aluminum  near  the  endurance 
limit),  the  maximum  effect  is  a  factor  3.4;  the  actual  nonlinear  life¬ 
time  would  be  (1/3.4)  «  0.3  times  the  linear  estimate  of  lifetime.  We 
note  that  these  effects  are  small  when  expressed  in  terms  of  the  es¬ 
timated  response;  the  maximum  effect  of  nonlinearity  on  lifetime  corres¬ 
ponds  to  an  Increase  of  8^  (0.7  db)  in  the  rms  response. 

The  three  curves  labelled  CASE  A  in  Pig.  2  show  the  effect  of  non¬ 
linearity  on  damage  rate  when  the  membrane  strain  (the  nonlinear  term 
in  Eq.  (19))  is  neglected  and  the  fatiguing  strain  is  linearly  propor¬ 
tional  to  displacement.  Because  nonlinearity  diminishes  the  probability 
of  large  displacements,  the  damage  rate  is  uniformly  diminished  in  Case  A. 
Asymptotic  expressions  for  the  ratios  of  deimage  rates  for  Cases  A  and  B, 
valid  for  extreme  nonlinearity,  have  been  computed  and  published  (ref.  I8). 
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AVERAGE  DAMAGE  PER  CYCLE  NORMALIZED  TO  THE  LINEAR  ESTIMATE 


0.5  1.0  10  100 

LINEARITY  PARAMETER,  Y 


FIG.2  AVERAGE  RATE  OF  ACCUMULATION  OF 
DAMAGE  (PER  CYCLE)  IN  RANDOM 
RESPONSE  OF  A  HARD-SPRING  RESONATOR 
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The  curves  In  Fig.  2  only  show  changes  from  the  linear  prediction 
of  the  damage  rate.  The  linear  prediction  of  damage  rate  was  discussed 
In  subsection  A  above.  It  Is  encompassed  In  the  "equivalent  strain 
level"  [Eq.  (4)],  which  Is  defined  as  the  steady  response  level  that 
has  the  same  lifetime  as  the  specified  rms  strain  response  to  random 
excitation.  The  difference  (in  decibels)  between  the  equivalent  strain 
level  and  the  rms  strain  level  Is  plotted  In  Fig.  3  against  a,  the  slope 
of  the  log  e  -  log  N  curve. 

2.  Range  of  Damaging  Strains 

We  consider  here  ways  of  characterizing  the  ranges  of  values  of 
strain  peaks  which  are  responsible  for  the  greatest  part  of  the  fatigue 
damage. 

In  the  linear  case  this  Information  Is  readily  shown  by  plotting 
the  cumulative  distribution  function  of  fatigue  damage.  These  curves 
show  as  ordinate  the  percent  of  total  damage  which  Is  attributable  to 
response  strain  peaks  smaller  than  the  abscissa  value.  Such  curves  are 
shown  In  Fig.  k,  where  the  abscissa  Is  the  strain  level  In  db  re  the  rms 
value  of  strain.  For  consistency  with  the  nonlinear  results.  Fig.  4  Is 
based  on  assumptions  of 

-  Miner's  law  of  damage  accumulation  [Eq.  (5)] 

-  straight  fatigue  curve  In  a  log-log  plot  [Eq.  (3)] 

-  Rayleigh  distribution  of  strain  peaks  [Eqs.  (l)  and  (2)]. 

The  resulting  damage  distribution  function  Is  an  Incomplete  Gainma  func¬ 
tion,  for  which  tables  exist.* 

Also  plotted  on  Fig.  4  are  the  points  of  maximum  damage  density. 
These  are  the  points  at  which  the  damage  density  —  the  Integrand 

x°^p(x)  =  exp(-x^/2) 

in  the  Integral,  Eq.  (5)^  for  total  damage  —  has  Its  maximum  value. 
(Thus  It  Is  the  "mode"  of  the  damage  density  curve.)  Analytically,  the 
point  of  maximum  damage  (sometimes  called  "peak  damage  point")  Is 


*The  function  plotted  In  Fig.  4  Is,  In  the  notation  of  Pearson's  Tables 
(ref.  19), 
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FIG.3  LINEAR  PREDICTION  OF  EQUIVALENT 
STRAIN.  STEADY  RESPONSE  ATEQUIV- 
ALENT  STRAIN  LEVEL  YIELDS  SAME 
LIFETIME  AS  RANDOM  RESPONSE  AT 
RMS  STRAIN  LEVEL 
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LEVEL  OF  PEAK  OF  EQUIVALENT  STRAIN 
DB  RE  RMS  RANDOM  STRAIN 
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FIG.4  CUMULATIVE  DISTRIBUTION  CURVES  FOR  DAMAGE  (LINEAR  CASE) 


(21) 


XpD  =  {a+1)^^^ 

where  x  =  ^/o,  is  the  ratio  of  strain  to  the  rms  value  of  strain. 

This  point  is  seen  to  lie  very  close  to  the  median  of  damage,  i.e., 
the  50$^  point  in  Fig.  4. 

Comparable  detail  about  the  distribution  of  damage  in  the  non¬ 
linear  random  response  is  not  available.  However  some  comparisons 
between  linear  and  nonlinear  cases  can  be  made.  The  most  important 
range  of  values  of  response  strain  is  conveniently  characterized  by 
the  values  of  the  quartile  points  on  the  distribution  curves,  the 
points  which  have  ordinate  values  of  25^  and  75/^  on  Fig.  4.  We  shall 
use  here  the  difference  in  the  levels  of  the  quartile  strains  (that  is, 
20  times  the  logarithm  of  their  ratio)  as  a  measure  of  the  range  of 
damaging  strains.  It  was  found,  by  hand  calculations  from  the  results 
of  the  digital  computer  analysis,  that  this  range  of  damaging  strains 
is  essentially  constant,  as  shown  in  Table  I.  The  differences  shown 
(0.3  db  or  less  than  are  not  significant  in  view  of  the  possible 
errors  in  hand  manipulation  of  the  data. 

TABLE  I 


THE  INTERQUARTILE  RANGE  OF  DAMAGE,  IN  DECIBELS:  20  . 


Values  for  nonlinear  response 
parameter  y  [Eq.  (20)]  from  4 

were  computed  for 
to  32,  and  varied 

values  of  linearity 
no  more  than  +0.1  db 

slope  a  of 
fatigue  curve: 

4 

8 

16 

linear  range 

3.6 

2.8 

2.0 

nonlinear  range 

3.9 

2.9 

2.1 

As  a  measure  of  central  tendency  of  the  daunage  distribution 
curves  in  the  case  of  nonlinear  response,  we  shall  use  the  geometric 
mean  value  of  the  strains  corresponding  to  the  quartile  points  (25/^ 
and  75^)  on  the  curve.*  In  the  linear  case  the  geometric  mean  of  the 
quartile  points  is  very  close  (within  0.1  db  or  1^)  to  the  median 
(50^  point).  In  the  presence  of  response  nonlinearity,  the  geometric 
mean  value  is  slightly  larger  than  for  the  corresponding  linear  struc 
ture  (Fig.  5).  The  amount  of  change  correlates  quite  well  with  the 
change  in  damage  rate  shown  in  Fig.  2  (Case  B) . 


*This  somewhat  unusual  choice  was  dictated  by  the  availability  of  data. 
Values  of  the  median  of  the  nonlinear  damage  curve  were  not  recorded 
in  the  course  of  computations.  The  point  of  maximum  damage  density 
(the  mode)  suffers  from  the  disadvantage  of  being  sensitive  to  the  choice 
of  independent  variable.  However  some  values  are  given  in  Fig.  5  of 
reference  I6. 
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GEOMETRIC  MEAN  OF  QUARTILE  DAMAGE  STRA 


FIG.  5  GEOMETRIC  MEAN  OF  THE  STRAINS  AT 
THE  QUARTILE  POINTS  OF  DAMAGE 
DISTRIBUTION  CURVE,  AS  A  FUNCTION 
OF  NONLINEARITY 
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3.  Summary 


An  analysis  was  made  of  the  effects  of  dynamic  nonlinearity  upon 
fatigue  damage  from  random  excitation.  The  results  are  strictly  per¬ 
tinent  only  to  a  particular  model  of  the  structure:  A  beam  with  pinned 
ends,  vibrating  In  Its  fundamental  mode,  with  fatigue  damage  occurring 
on  the  beam's  surface  at  midspan.  However,  It  Is  expected  that  the 
present  results  would  be  qualitatively  pertinent  to  other  nonlinear 
structures  wherein  the  response  nonlinearity  arises  from  a  similar 
cause,  namely  a  symmetric  but  nonlinear  spring  force.  (Other  struc¬ 
tures  with  similar  dynaunlcal  equations  are  discussed  In  a  recent  review 
paper  (ref.  20).) 

It  was  found  that,  with  moderate  nonlinearity,  the  rate  of  accu¬ 
mulation  of  damage  Is  Increased  but  slightly  from  the  rate  one  would 
estimate  by  Ignoring  the  nonlinearity  (Fig.  2).  The  change  In  damage 
rate  depends  upon  the  extent  of  nonlinearity  and  the  slope  of  the 
fatigue  curve.  For  values  of  slope  from  4  to  l6,  the  change  In  damage 
rate  corresponds  to  a  change  In  excitation  not  exceeding  0.7  db  (8j6). 

Secondly  It  was  found  that  nonlinearity  causes  very  little  change, 
from  the  linear  estimate.  In  the  absolute  values  of  the  strains  that  do 
the  most  damage.  This  fact  was  shown  both  by  the  essential  constancy 
In  the  ratio  of  the  quartlle  strains  In  the  damage  distribution  (Table  I) 
and  by  the  small  changes  In  the  geometric  mean  of  the  quartlle  strains 
(Fig.  5)*  These  latter  changes  (small  Increases)  have  about  the  right 
magnitude  to  account  for  the  changes  In  damage  rate  shown  In  Fig.  2. 


C.  EXPERIMENTAL  EQUIVALENCE  SCHEME  FOR  SYSTEM 
WITH  RESPONSE  NONLINEARITY  AND  SEVERAL  MODES 

An  equivalence  scheme  using  experimentally-determined  engineering 
corrections  for  nonllnearltles  of  response,  and  for  the  presence  of 
more  than  one  mode,  has  been  outlined  In  a  number  of  reports  (ref.  4). 
Basically  this  procedure  Is  an  adaptation  and  modification  of  the  scheme 
outlined  In  subsection  A  above. 

Suppose  that  It  Is  required  that  a  structure  withstand  a  certain 
Jet  noise  environment  (specified  spectral  density).  A  structure  has 
been  designed  and  samples  of  It  built.  One  wishes  to  specify  response 
and  fatigue  tests  using  a  pure  tone  siren  and  a  computation  procedure 
that  will  determine  the  average  life  to  be  expected  of  the  structure 
In  the  given  environment.  The  determination  would  proceed  in  the 
following  manner.* 


*It  Is  not  entirely  evident  from  the  reports  (ref.  4)  that  the  test 
procedure  described  here  Is  that  which  was  actually  used  In  practice. 
However  we  have  attempted  here  to  be  faithful  to  the  principles  described 
In  the  references. 
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1.  Response  Tests  and  Linear  Prediction  of  Random  Response 


By  a  combination  of  prior  experience,  numerous  measurements 
(including,  perhaps,  strain  gauges  in  many  places,  strain-sensitive 
lacquer,  and  qualitative  fatigue  tests  using  the  siren),  and  good 
luck,  a  single  location  and  orientation  for  a  reference  strain  gauge 
is  chosen;  its  output  represents  the  nominal  response  strain  used  in 
calculations.  The  response  of  the  structure  to  pure  tone  (siren) 
excitation  is  investigated  as  a  function  of  frequency,  to  detect  the 
various  modes  of  response,  and  as  a  function  of  sound  pressure  level 
to  detect  nonlinearity  in  the  various  modes.  The  latter  tests  are 
made  with  sound  frequencies  equal  to  the  resonance  frequencies  of  the 
various  modes.  Generally  it  is  found  that  motion  in  one  of  the  modes 
predominates . * 

A  linear  prediction  of  the  rms  nominal  strain  In  the  random  environ¬ 
ment  is  made  by  the  techniques  discussed  in  subsection  A  Eq. (17),  based 
on  the  results  of  the  response  measurements.  This  prediction  is  made  for 
each  important  mode. 

2.  Correction  for  Extra  Modes 


When  more  than  one  mode  is  apparent  in  the  response,  the  predicted 
mean-square  responses  for  the  various  modes  are  simply  added  to  yield  a 
prediction  of  the  total  mean-square  response  in  the  random  environment. 

3.  Correction  for  Response  Nonlinearity 

The  linear  prediction  equation,  Eq.  {l7)t  for  the  rms  of  nominal 
strain  in  random  response  depends  upon  the  ratio  of  nominal  strain  to 
sound  pressure  which  is  found  by  pure  tone  response  measurements  at 
resonance.  When  the  response  of  a  mode  is  nonlinear,  that  ratio  varies 
with  strain  level. 

In  the  present  equivalence  scheme,  compensation  is  made  for  such 
nonlinearity  by  requiring  that  the  strain-pressure  ratio  be  evaluated 
at  the  "peak-damage  strain",  l.e.,  the  level  of  nominal  strain  which 
lies  at  the  maximum  of  the  daunage  density  curve.  (See  subsection  B-2, 
Eq.  (21),  and  Fig.  4,  above.)  Of  course  thlb  level  cannot  be  predicted 
exactly  at  this  stage,  since  the  ratio  of  peak-damage  strain  to  rms 
strain  depends  on  lifetime  and  the  lifetime  depends  on  the  strain  con¬ 
centration  factor  which  has  not  been  determined.  However  a  fairly  good 
estimate  can  usually  be  made. 


*In  the  present  context,  "mode"  refers  to  a  single,  narrow  peak  in  the 
response  curve. 
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After  a  strain  concentration  factor  has  been  determined  by  pure- 
tone  fatigue  testing  (subsection  h,  below),  a  corrected  estimate  of 
the  value  of  the  peak-damage  strain  is  made  and  the  prediction  of  rms 
random  response  is  corrected  according  to  the  change  in  nonlinearity 
of  pure-tone  response. 

The  order  in  which  corrections  should  be  made  for  extra  modes  and 
for  nonlinearity,  and  the  manner  in  which  these  corrections  may  inter¬ 
act,  are  not  entirely  clear.  When  one  mode  is  dominant,  these  problems 
are  probably  of  little  importance.  The  correction  for  extra  modes 
should  then  be  a  small  correction. 

4.  Fatigue  Tests  and  Prediction  of  Life  in  Random  Environment 

From  constant-amplitude  mechanical  fatigue  tests  on  specimens  of 
the  material,  designed  with  no  strain  concentrations,  a  fatigue  curve  is 
determined.  These  data  may  come  from  the  available  literature. 

A  constant-amplitude  fatigue  test  is  made  on  the  structure,  excited 
by  sound  from  the  siren,  at  a  (peak)  nominal  strain  equal  to  the  expected 
peak-damage  nominal  strain.  The  lifetime  is  measured.  The  difference 
between  the  experimental  nominal  strain  and  the  strain  on  the  fatigue 
curve  which  corresponds  to  the  seime  lifetime  represents  the  strain  con¬ 
centration  factor.  This  factor  is  applied  uniformly  to  the  old  fatigue 
curve  to  generate  a  new  fatigue  curve  relating  nominal  strain  and  life¬ 
time  . 


On  first  thought,  it  might  seem  desirable  to  try  to  generate  a 
fatigue  curve  for  nominal  strain  by  repeated  fatigue  tests  on  nominally 
Identical  structures.  However,  the  structural  variability  which  must  be 
expected  between  specimens  and  the  relative  awkwardness  of  the  leviathan 
which  is  a  siren  "facility"  militate  against  such  a  procedure. 

From  the  constant-amplitude  fatigue  curve  for  nominal  strain  which 
has  now  been  constructed,  the  lifetime  for  random  response  with  the  pre¬ 
dicted  rms  nominal  strain  is  computed  by  the  procedures  of  subsection 
A-1  above  [Eq.  (5)]^  assuming  a  Rayleigh  distribution  of  strain  peaks 
and  ignoring  any  fatigue  interaction  effects. 

5.  Comments 


The  equivalence  scheme  described  in  this  subsection  has,  through 
its  use  in  the  aircraft  Industry,  been  proved  a  practical  engineering 
procedure  although  some  questions  have  been  raised  over  its  accuracy. 
It  is  evident  that  no  attempt  is  made  to  compensate  for  any  fatigue 
interactions  between  response  peaks  of  different  amplitude. 

The  correction  for  the  presence  of  more  than  one  response  mode  is 
very  much  an  hoc  procedure  with  little  analytical  basis.  When  the 
several  resonant  response  modes  are  well  separated,  the  procedure  of 
summing  mean-square  responses  in  the  separate  modes  yields  correct 
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predictions  of  the  total  mean-square  response.  However,  there  is  no 
available  information  to  Indicate,  one  way  or  the  other,  whether  the 
relationship  between  fatigue  life  and  total  mean-square  response  is 
the  same  for  one  and  for  several  resonant  modes.  Use  of  the  correc¬ 
tion  is  at  least  a  more  conservative  procedure  than  the  Ignoratlon  of 
extra  modes. 

The  correction  for  nonlinearity  can  be  based  on  two  hypotheses 
concerning  the  random  nonlinear  response  mechanism: 

(1)  Insofar  as  significant  fatigue  damage  is  concerned,  the  non¬ 
linear  distribution  of  strain  peaks  can  be  represented  as  a  Rayleigh 
distribution  with  a  new  value  of  rms  strain.  Thus,  If  the  peak  den¬ 
sity  curve  Is  plotted  on  a  logarithmic  scale  (e.g.,  linear  In  decibels), 
the  effect  of  nonlinearity  Is  merely  a  displacement  of  the  curve  without 
distortion. 

(2)  The  extent  of  this  displacement  Is  equal  to  the  deviation  from 
linearity  of  the  strain-pressure  ratio,  measured  at  pure-tone  resonance 
at  the  "peaik- damage"  strain. 

The  exact  nonlinear  analysis  given  In  the  previous  subsection  can 
be  used  to  check  these  hypotheses,  for  the  case  of  one  special  type  of 
nonlinearity.  The  validity  of  the  first  hypothesis  is  fairly  well 
borne  out  by  the  observed  near-constancy  of  the  Inter-quartlle  range 
of  daunage  [subsection  A-2  and  Table  I].  The  validity  of  the  second 
hypothesis  has  not  been  checked  in  detail.  A  few  computations  have  been 
compared  with  the  results  given  In  Pig.  5-  They  show  that  the  second 
hypothesis  gives  the  correct  direction  for  the  displacement  of  the  damage 
curve,  but  underestimates  the  magnitude  of  the  displacement  (in  logarithmic 
measure,  l.e.,  in  decibels)  by  a  factor  of  about  1/2.  [The  calculations 
are  made  from  Eqs.  (7)  and  (ll)  of  reference  I6.  ] 

However,  it  must  be  noted  that  the  hard- spring  model  of  a  nonlinear 
structure  used  in  these  calculations  Is  a  very  idealized  and  specialized 
case  of  nonlinearity..  Other  types  of  response  nonlinearity  may  be  pre¬ 
sented  and  may  dominate;  nonlinearity  of  damping  is  particularly  notable 
in  experimental  measurements  (ref.  21).  There  Is  no  analytical  basis 
for  verifying  these  hypotheses  in  the  case  of  nonlinear  damping. 


D.  EXPERIMENTAL  EQUIVALENCE  SCHEME  FOR  SYSTEM  WITH 
RESPONSE  NONLINEARITY  AND  FATIGUE  INTERACTIONS 

We  outline  In  this  section  a  new  slren-Jet  fatigue  equivalence 
scheme  which  does  not  depend  upon  calculations  based  on  the  absence  of 
fatigue  Interactions,  l.e.,  on  Miner's  rule  [Eq.  (5)].  Moreover,  the 
effects  of  nonlinearity  of  response  are  Included. 
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The  proposed  procedure  is.  In  essence,  to  use  an  experimentally- 
determined  programmed  fatigue  teat  In  which  the  structure  Is  driven 
at  resonance  by  a  sound  pressure  whose  amplitude  varies  In  time  accord¬ 
ing  to  a  predetermined  program.  The  expected  total  test  time  Is  divided 
Into  a  number  of  equal,  nominally  Identical  blocks.  Within  each  block 
the  excitation  and  response  Increase  to  a  maximum  value  and  then  decrease 
to  zero  (or  a  close  approximation  thereto).  The  variations  may  be  con¬ 
tinuous  or  In  steps  chosen  to  approximate  the  proper  continuous  curve. 

The  program  of  pressure  Is  to  be  chosen  so  that  Its  probability 
density  reproduces  the  probability  density  of  strain  In  random  response 
(except  for  a  uniform  scale  change)  when  the  response  Is  linear.  Thus 
If  the  overall  excitation  Is  small  so  that  response  Is  linear,  the  pres¬ 
sure  progreim  Is  chosen  so  that  the  long-time  (one  block)  probability 
density  of  strain  Is  Identical  with  the  density  of  strain  under  random 
excitation  by  some  one  value  of  rms  random  sound  pressure. 

Experimentally,  the  program  can  be  designed  as  follows.  The  struc¬ 
ture  Is  exposed  to  a  low  level  of  random  noise  while  It  Is  mounted  In  an 
environment  essentially  Identical  to  the  Jet-noise  environment.  Measure¬ 
ments  are  made  of  the  rms  sound  pressure,  the  rms  strain  response,  and 
the  density  of  peaks  In  that  response. 

Then  the  structure  Is  mounted  In  position  for  the  Intense  pure-tone 
(l.e.,  siren)  fatigue  test.  It  Is  excited  by  a  low-level  pure  tone. 

The  sound  pressure  Is  adjusted  until  the  rms  strain  In  this  test  equals 
the  rms  strain  In  the  previous  random  test.  The  rms  value  of  the  adjusted 
sound  pressure  Is  defined  as  the  equivalent  of  the  rms  sound  pressure  In 
the  random  test. 

Now,  with  linear  response,  the  pure-tone  response  Is  proportional 
to  the  pure-tone  sound  pressure.  Therefore  the  desired  program  of  sound 
pressures  Is  one  having  the  same  peak  statistics  (relative  to  the  equiva¬ 
lent  rms  pressure)  as  the  peak  statistics  of  the  random  strain  (relative 
to  the  rms  random  strain).  Analytical  and  practical  details  are  discussed 
In  section  V  below. 

By  following  these  steps  we  have  determined  a  program  of  low-level 
pure-tone  sound  pressures  which.  In  the  siren  environment.  Is  equivalent 
to  a  measured  rms  random  sound  pressure  In  the  Jet  environment.  We  postu¬ 
late  that  the  equivalence  of  responses  will  be  maintained  If  the  pure-tone 
and  random  sound  pressures  are  all  scaled  up  by  the  same  factor  (l.e.,  all 
sound  pressure  levels  Increased  by  the  same  number  of  decibels).  We  fur¬ 
ther  postulate  that  the  fatigue  lifetimes  In  siren  and  Jet  environments 
will  be  equal.  If  the  scaDe  factor  brings  the  response  Into  the  fatigue 
region. 

One  of  the  major  purposes  of  this  study  Is  to  check  these  postulates 
experimentally.  There  are  two  underlying  aspects  to  them  that  can  be 
discussed  before  the  data  are  presented.  These  are  (l)  the  Influence  of 
any  fatigue  damage  Interactions  between  strain  peaks  of  different  ampli¬ 
tudes,  and  (2)  the  Influence  of  response  nonlinearities.  Those  discus¬ 
sions  follow  the  next  paragraph. 


ASD  TR  61-639 


28 


The  preceding  exposition  of  procedure  for  determining  the  pressure 
program  was  made  simple  by  ignoring  any  practical  problems.  Preliminary 
testing  to  help  choose  a  position  for  the  strain  gauge  is  as  important 
as  in  the  preceding  scheme  (subsection  II-C).  The  choice  between  meas¬ 
uring  pressures  with  the  panel's  motion  blocked  or  with  the  panel  moving 
must  be  made  (subsection  II-A).  Moreover  an  experimenter  would  certainly 
measure  the  ratio  of  pure-tone  pressure  to  strain  at  more  than  one  ampli¬ 
tude.  Further  belaboring  of  such  points  is  unnecessary  for  experimentally 
inclined  readers  and  unrewarding  to  the  theoretician. 

1.  Fatigue  Interactions 

Fatigue  experiments  made  in  the  laboratory  with  strain  cycles  of 
varying  amplitude  have  long  been  notorious  for  their  disagreement  with 
predictions  of  lifetime  based  on  so-called  "linear  damage  laws",  l.e., 
theories  that  neglect  fatigue  Interactions  between  cycles  of  different 
amplitudes.  Freudenthal ' s  results,  for  example,  have  demonstrated  the 
importance  of  such  effects  in  his  tests  on  aluminum  alloys  2024  and  7075 
and  on  SAE  4340  steel  alloy  where  he  generally  used  a  random  response 
with  an  exponential  density  of  stress  peaks  (refs.  8,  17^  and  22-25). 

The  various  formalisms  which  have  been  proposed  to  match  the  experi¬ 
mental  data  —  fictitious  S-N  curves,  etc.  —  do  not  take  a  form  which  is 
generally  useful  for  predicting  fatigue  life,  for  a  variety  of  reasons. 
Some  Involve  parameters  not  experimentally  verifiable.  In  other  cases 
the  dependence  of  the  parameters  upon  the  distribution  of  the  random 
loading  is  unknown.  Such  a  readily  adaptable  procedure  as  the  fictitious 
S-N  curve,  with  a  constant  slope  in  a  log-log  plot,  may  yield  a  conserva¬ 
tive  estimate  of  fatigue  lifetime,  but  it  may  also  be  so  conservative  as 
to  constitute  an  unwanted  design  penalty.  Formalisms  dependent  on  a 
value  called  "maximum  stress"  lead  to  special  frustration  when  only  sta¬ 
tistical  probabilities  for  the  occurrence  of  given  response  levels  are 
known . 

In  the  present  slren-Jet  fatigue  equivalence  scheme,  we  have  pro¬ 
posed  to  use  the  method  of  "programmed  loading"  to  Incorporate  the 
effects  of  fatigue  interactions  in  the  siren  tests. 

Unfortunately  the  literature  does  not  present  a  single  picture  of 
the  success  of  such  programmed  loading  in  reproducing  the  fatigue  effects 
of  random  loading.  Where  some  tests  of  programmed  loading  yield  shorter 
lifetimes  than  random  loading  (ref.  22,  Fig.  2),  others  show  the  opposite 
trend  (ref.  26,  Fig.  l4).  These  differences  may  very  well  be  related  to 
differences  in  the  distributions  of  the  random  strains  and  to  differences 
in  the  short-time  correlation  between  strains  of  different  levels  in  the 
random  signal. 
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2.  Response  Nonlinearity 


We  now  consider  the  effect  of  response  nonlinearity  on  the  postulate 
that  the  equivalence  between  progranmed  response  and  random  response  is 
invariant  under  a  uniform  scaling  of  the  excitations.  First  consider  the 
case  of  a  structure  with  nonlinear  stiffness,  for  which  exact  statistical 
solutions  are  available. 

The  analysis  for  random  response  in  one  mode  with  a  nonlinear  stiff¬ 
ness  was  summarized  in  subsection  II-B  above.  The  following  points  are 
pertinent.  The  probability  density  of  the  velocity  (time  rate  of  change 
of  displacement)  is  Gaussian  and  its  rms  value  is  linear  in  the  strength 
of  excitation  (ref.  27).  Therefore  the  density  of  velocity  peaks  is  a 
Rayleigh  distribution  at  all  levels  of  excitation.  The  density  of  dis¬ 
placement  peaks  is  not  Rayleigh.  However,  if  the  relation  between  peak 
displacement  and  peak  velocity  which  is  valid  for  pure-tone  steady-state 
vibration  at  resonance  is  used  to  transform  the  (Rayleigh)  density  of 
velocity  peaks,  then  the  correct  density  of  displacement  is  obtained 
(ref.  20). 

Now,  the  peak  velocity  of  nonlinear  response  to  pure-tone  steady- 
state  excitation  at  resonance  is  linear  in  the  strength  of  excitation 
(ref.  21).  If,  therefore,  the  nonlinear  structure  is  excited  to  quasl- 
steadystate  resonance  by  a  pure-tone  excitation  whose  strength  is  varied 
so  as  to  result  in  a  Rayleigh  distribution  of  peak  values,  then  the 
long-time  density  of  its  displacement  peaks  will  have  the  same  distribu¬ 
tion  as  in  the  case  of  random  excitation.  Since  the  peak  strain  is 
uniquely  determined  by  peak  displacement,  the  peak  strains  must  also 
have  the  same  statistics  with  random  and  programmed  excitations.  This 
completes  the  proof  of  the  equivalence  scheme  for  structures  with  non¬ 
linear  stiffness. 

Other  types  of  nonlinearity,  for  example  in  damping,  are  significant 
in  many  structures,  as  we  remarked  in  subsection  II-C.  There  is  no  anal¬ 
ysis  by  which  the  effects  of  such  nonlinearity  on  the  response  statistics 
can  be  assessed.  However  the  random  response  of  a  highly  resonant  struc¬ 
ture,  on  the  average,  probably  resembles  closely  a  sequence  of  natural 
vibration  cycles  whose  peak  amplitude  varies  slowly  in  time.  (See  the 
discussion  in  ref.  20.)  Retention  of  the  present  equivalence  scheme  under 
these  circumstances  is  intuitively  attractive  until  some  pertinent  analyt¬ 
ical  or  experimental  evidence  to  the  contrary  is  found. 
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SECTION  III 


FATIGUE  TESTS  OP  RESONANT  STRUCTURE  UNDER 
RANDOM  AND  CONSTANT- AMPLITUDE  EXCITATIONS 


A  series  of  fatigue  tests  were  made  with  a  simple  resonant 
structure  excited  mechanically  by  random  forces  and  by  sinusoidal 
forces  of  constant-amplitude.  The  basic  design  of  these  tests  is 
very  similar  to  those  used  by  Head  and  Hooke  (ref.  3)  and  Prallch 
(ref.  5). 


The  majority  of  these  tests  were  made  on  beam  specimens  of  plain 
2024  T3  aluminum  alloy,  sheet,  excited  to  bending  vibrations.  The 
results  of  the  constant-amplitude  tests,  with  failure  occuring  at  a 
transverse  hole,  are  compared  with  constant-aimplitude  tests  reported 
in  the  literature,  to  determine  a  value  of  strain  concentration  due 
to  the  hole.  The  constant-5implltude  results  are  also  combined  ac¬ 
cording  to  Miner's  rule  [Eqs.  (l)  and  (5),  subsection  II-A]  to  yield 
predictions  of  fatigue  life  with  random  excitation  which  are  compared 
with  the  results  of  the  random  experiments. 

Other  tests  were  made  with  random  excitation  of  similar  specimens 
of  Alclad  2024  aluminum  alloy  to  determine  a  relation  between  fatigue 
lifetime  and  rms  response.  This  relation  will  later  be  compared  with 
results  of  excitation  by  Jet  noise  of  the  larger  structural  panels, 
which  were  made  of  this  material. 


A.  APPARATUS  AND  TEST  PROCEDURES 

The  test  specimens  are  simple  rectangular  beaims  cut  from  commer¬ 
cial  sheet  stock  (O.O63  in.  thick)  to  the  dimensions  given  in  Pig.  6. 
For  the  tests,  a  specimen  is  clamped  across  the  center  between  two 
l/4-in.  bars  which  form  the  support  by  which  it  is  forced.  The  two 
middle  holes  in  Pig.  6  are  positioning  holes  covered  by  this  support. 
Thus  the  specimen  and  the  supporting  clamp  form  a  double  cantilever 
beam  system.  When  the  clamp  is  forced  at  a  suitably  chosen  frequency 
the  two  ends  of  the  specimen  vibrate  with  large  amplitudes  while  the 
clamp  Itself  moves  very  little.  The  response  is  principally  in  the 
fundamental  cantilever  mode,  with  a  slight  amount  of  rigid  body  dis¬ 
placement  superposed. 

The  two  1/8-ln .  holes  drilled  through  the  specimen  at  distances 
of  0.5  in.  from  the  transverse  center  line  (Pig.  6),  provide  con¬ 
trolled  strain  concentrations.  They  are  necessaiy  to  prevent  un¬ 
observable  fatigue  failures  at  the  edge  of  the  clamp;  moreover,  fail¬ 
ure  in  bending  at  a  hole  is  a  common  occurrence  in  aircraft  structures 
(e.g.,  at  rivet  holes). 

The  specimens  were  made  with  normal  care  for  dimensions,  etc., 
in  our  machine  shop.  Batches  were  milled  and  drilled  together;  no 
variations  between  batches  were  noted  in  the  fatigue  tests.  The 
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FIG.6  CANTILEVER  TEST  SPECIMEN  FOR  RANDOM  AND 
CONSTANT-AMPLITUDE  FATIGUE  STUDIES 
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long  side  of  the  beams  was  parallel  to  the  direction  in  which  the 
plate  stock  was  rolled.  The  strain  concentrating  holes  were  drilled 
to  size  in  stages,  with  care  taken  to  avoid  burrs,  however  no  attempt 
was  made  to  polish  or  otherwise  finish  them. 

Instrumentation  and  test  procedures  are  Illustrated  in  Fig.  ?• 

The  central  supporting  clamp  of  the  test  specimen  is  screwed  securely 
to  the  voice  coll  of  a  Goodmans  Type  390A  shaker,  driven  electrically 
by  the  chosen  signal.  A  nominal  strain  response  Is  measured  by  a 
strain  gauge  In  the  transverse  plane  containing  one  of  the  strain  con¬ 
centrating  holes. 

The  response  level  Is  also  determined  by  measuring  the  displace¬ 
ment  of  the  end  of  the  beam.  The  value  of  strain  calculated  from  the 
displacement,  using  the  theoretical  shape  for  the  lowest  mode  of  a 
cantilever,  was  found  to  be  In  good  agreement  with  that  measured  by 
the  strain  gauge,  after  a  correction  was  made  for  the  amplitude  of 
clamp  displacement.  Measurement  of  displacement  was  found  to  be  es¬ 
sential  for  monitoring  strain  level  during  fatigue  tests  because  the 
strain  gauges  usually  failed  after  a  few  minutes  at  the  higher  test 
levels . 

In  fatigue  tests  with  random  excitation,  electrical  noise  from  a 
random  noise  generator  Is  passed  through  a  band-pass  filter  tuned  to 
a  broad  frequency  region  centered  on  the  fundamental  resonance  fre¬ 
quency  of  the  beams.  This  procedure  limits  low-frequency  signals 
that  might  cause  too  large  a  displacement  of  the  voice  coll  and  limits 
the  amplitude  of  response  In  the  higher  cantilever  modes.  A  SOO-watt 
power  amplifier  drives  the  shaker. 

Basic  measurements  of  strain  response  are  made  with  a  Ballantlne 
Model  320  true  rms  voltmeter,  A  long  averaging  time  (8  sec  time 
constant)  Is  required  to  obtain  reasonably  precise  measurements  be¬ 
cause  of  the  narrow  bandwidth  of  the  beam's  resonance  (Q  »  45,  c/c^,  * 
0.01).  So  that  the  beam  should  not  suffer  significant  fatigue  damage 
while  the  response  level  Is  being  set  and  measured.  It  Is  necessary 
to  make  the  measurement  at  a  slightly  lower  level  ajid  then  Increase 
the  excitation  by  a  known  amount  with  an  accurate  attenuator.  Tape 
recordings  of  the  random  response  strains  are  made,  as  long  as  the 
strain  gauge  lasts,  for  subsequent  statistical  analysis. 


B.  RESULTS  AND  ANALYSIS 
1 .  Constant-Amplitude  Tests 

The  results  of  these  fatigue  measurements  on  the  beams  of  plain 
2024  aluminum  alloy  are  summarized  In  Fig.  8.  Three  samples  were 
fatigued  at  each  of  four  levels  with  constant-amplitude  vibrations. 
Through  these  points  Is  drawn  the  smooth  curve  B. 
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FIG.  7  INSTRUMENTATION  FOR  CANTILEVER  BEAM  FATIGUE  TEST 
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FIG.  8  FATIGUE  CURVES  FOR  2024  ALUMINUM.  RESULTS  OF  CONSTANT-AMPLITUDE 
AND  RANDOM  EXPERIMENTS  ON  RESONANT  CANTILEVER  BEAM,  AND 
PREDICTIONS  THEREFROM. 


The  fatigue  curve  A  Is  determined  by  the  experimental  data  ob¬ 
tained  by  Freudenthal  (ref.  8)  and  D' Amato  (ref.  7)  for  test  specimens 
\vlthout  marked  concentration  factors.  Freudenthal '  s  data  were  acquired 
by  the  rotatlng-beam  method.  His  values  of  stress  were  converted  to 
strain  by  the  Young's  modulus:  1.06  x  107  psl.  D' Amato's  results, 
taken  In  tensile  tests,  were  used  to  extend  the  curve  beyond  the  yield 
point,  for  reasons  explained  In  subsection  II-A-2,  above. 

The  difference  between  the  curves  Is  an  experimental  strain  con¬ 
centration  factor  for  fatigue.  In  logarithmic  measure.  The  values 
vary  from  4.1 'db  at  the  lowest  test  level  to  3-0  db  at  the  highest 
level.  The  classical  elastic  stress  concentration  factor  for  a  hole 
In  a  very  wide  plate  In  bending  (Poisson's  ratio  of  1/3)  Is  1.9,  cor¬ 
responding  to  5-5  db  (r3f.  28).  However,  Peterson  (ref.  29)  has  pro¬ 
posed  a  formula  for  estimating  the  effective  strain  concentration 
factor  for  fatigue: 

K.  +  a/r 

K„  =  - 

^  1  +  a/r 

where  K,  Is  the  elastic  concentration  factor,  a  Is  a  material  constant, 
and  r  Is  the  radius  of  the  hole.  For  2024  aluminum  alloy,  Peterson 
estimates  that  a  Is  0.05  In.  Therefore,  with  a  hole  radius  of  0.0625 
In.,  Peterson's  formula  yields  a  fatigue  concentration  factor  of  1.4, 
or  3-0  db  In  logarithmic  measure.  This  prediction  compares  quite  well 
with  the  experimental  values,  especially  at  the  higher  levels.  (The 
precision  of  measurement  of  strain  Is  estimated  to  be  +0.2  db.). 

It  Is  Interesting  to  note  that  the  fatigue  curve  obtained  by 
these  measurements  Is  parallel  to  and  approximately  0.5  db  lower  than 
the  results  of  earlier  tests  made  In  these  laboratories  on  the  fatigue 
life  of  a  cantilever  bar  held  by  a  line  of  Idealized  rivets  In  l/Q-ln. 
holes  (ref.  30) • 

2.  Random  Excitation  Tests 


Three  cantilever  test  specimens  were  fatigued  at  each  of  three 
rms  strain  levels;  the  results  of  these  tests  are  the  data  points 
labelled  E  In  Fig.  8.  The  probability  density  of  a  15  min.  sample 
of  the  nominal  strain  (about  1.6  x  105  cycles)  Is  plotted  In  Fig.  9. 
The  agreement  between  measurements  and  the  theoretical  density  of 
normal  (Gaussian)  noise  Is  quite  good. 

Fatigue  lifetime  as  a  function  of  the  rms  strain  was  computed 
from  the  constant-amplitude  fatigue  curves  In  Fig.  8  by  the  simple 
theory  of  subsection  II-A  above  LEqs.  (l)  and  (5)]^  assuming  a 
Rayleigh  distribution  of  the  strain  peaks  and  Ignoring  fatigue  Inter¬ 
actions.  The  resulting  random  fatigue  curve  for  plain  2024  aluminum 
alloy  without  strain  concentrations  Is  labelled  C.  The  random  fatigue 
curve  computed  from  the  experimental  constant-amplitude  data  Is 
labelled  D  In  Fig.  8. 
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PROBABILITY  DENSITY  OF  STRAIN  SIGNAL 


FIG. 9  PROBABILITY  DENSITY  OF  STRAIN  IN  CANTI¬ 
LEVER  BEAM  RESPONSE  TO  RANDOM  EXCITATION 
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It  will  be  noted  that  the  experimental  lifetimes  with  raindom 
excitation  are  shorter  than  the  prediction  which  ignores  fatigue 
interactions;  the  differences  are  greatest  at  the  longest  lifetimes. 
Freudenthal' s  detected  qualitatively  similar  differences  in  his  tests, 
made  with  the  same  material  but  with  no  strain  concentration  and  with 
a  response  which  differed  in  both  its  long-  amd  short-time  statistics. 
His  results  could  be  represented  fairly  well  by  a  simple  analytic 
expression  and  by  a  smooth  curve  showing  the  difference  between  ex¬ 
perimental  and  predicted  lifetimes  as  a  function  of  predicted  life¬ 
time  [ref.  17,  Eq.  (6.10)  and  Fig.  I6].  If  the  present  predicted  life 
times  (curve  D)  are  "corrected"  by  the  factor  given  in  Freudenthal ' s 
smooth  curve,  the  curve  labelled  F  in  Fig.  8  is  obtained.  The  present 
ejqperlmental  results  lie  between  the  two  curves  D  and  F. 

There  is  some  opportunity  for  uncertainty  inherent  in  the  detec¬ 
tion  of  a  fatigue  crack.  In  constant-amplitude  tests,  the  change  in 
resonance  frequency  and  the  consequent  drop  in  amplitude  if  frequency 
is  not  readjusted  are  very  sensitive  indicators  of  a  crack's  existence 
They  are  supplemented  by  visual  inspection  under  stroboscopic  light. 
However,  neither  of  these  techniques  is  adequate  when  the  excitation 
is  random;  then  the  excitation  must  be  turned  off  and  the  stationary 
specimen  Inspected  visually.  Thus,  in  the  random  fatigue  test,  the 
minimum  uncertainty  in  life  measurement  is  set  by  the  time  interval 
between  inspections. 

However,  the  rate  of  growth  of  crack  length,  measured  in  real 
time,  is  quite  slow  in  random  tests.  This  reflects  the  fact  that 
most  response  cycles  are  too  small  to  do  much  damage.  Some  measure¬ 
ments  made  on  one  specimen  are  shown  in  Fig,  10.  The  random  excita¬ 
tion  was  turned  off  from  time  to  time  and  the  specimen's  surface 
examined,  with  a  good  light  and  a  five-power  magnifying  glass.  The 
lengths  of  any  cracks  were  measured  by  comparison  under  the  glass 
with  a  rule  graduated  to  0.020  in.  In  Fig.  10,  these  data  are  plotted 
for  each  of  the  four  fatiguing  points  (the  four  edges  of  the  strain 
concentration  holes  which  lie  on  transverse  diameters). 

The  slow  rate  of  crack  growth  in  random  response,  possible  vari¬ 
ations  in  this  rate  with  level  of  response,  the  impossibility  of  con¬ 
tinuous  inspection,  and  the  dependence  of  the  visibility  of  a  given 
crack  upon  lighting  conditions  -  these  factors  can  create  quite  sig¬ 
nificant  errors  in  lifetime.  The  various  factors  may  bias  the  data 
in  different  directions,  favoring  early  detection  of  cracks  in  either 
the  long  or  the  short  tests.  They  may  Introduce  differences  between 
laboratory  and  field  data,  or  between  field  measurements  at  different 
times.  As  an  example,  in  one  field  test  reported  below,  two  speci¬ 
mens  were  exposed  to  the  same  noise  field;  the  test  was  concluded 
when  a  crack  was  apparent  in  one  specimen  and  not  the  other.  On 
laboratory  Inspection,  out  of  the  bright  sun,  noise,  and  other  dis¬ 
tractions,  one  could  readily  see  the  "missing"  crack  with  the  unaided 
eye. 
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GROWTH  OF  CRACK  LENGTH  DURING  STEADY  RANDOM 
EXCITATION  OF  RESONANT  BEAM 
(RESONANCE  FREQUENCY,  185  CPS) 


There  is,  therefore,  a  very  real  question  about  the  precision  of 
the  recorded  lifetimes,  and  not  much  to  guide  one's  Judgement.  How¬ 
ever,  the  trend  of  the  curves  in  Fig.  10  does  not  much  reassure  a 
natural  skeptic  of  the  reliability  of  a  50^  difference  in  life  re¬ 
corded  in  a  field  test. 


C.  DATA  FOR  ALCLAD  2024 

A  more  limited  investigation  of  the  same  kind  was  made  with 
cantilever  beam  specimens  of  Alclad  2024  aluminum  alloy.  The  design 
of  specimens  and  apparatus  was  the  same  as  in  the  earlier  tests  with 
plain  2024. 

The  purpose  of  these  experiments  is  to  determine  the  slope  of 
the  random  fatigue  curve  on  small  specimens  and  to  compare  this  slope 
with  the  slope  determined  in  the  Jet-noise  fatigue  tests.  Satisfactory 
correlation  of  results  will  indicate  that  the  small  specimen  tests 
could  have  been  used  for  prediction . 

Three  specimens  were  excited  by  steady  random  forces  at  each  of 
two  ms  strain  levels.  The  data  and  the  computed  average  slope  are 
given  In  Table  II. 


TABLE  II 

FATIGUE  DATA  FOR  CANTILEVER  BEAM 
SPECIMENS  OF  ALCLAD  2024 


Nominal 
Strain  Level 

(db  re  10"^) 

Specimen 

Number 

Fatigue 
Lifetime 
(10  log  N) 

56.0  db 

1 

59.2 

2 

59.2 

3 

58.9 

geometric  mean 

life: 

59.1 

61.0  db 

4 

49.4 

5 

49.8 

6- 

48.6 

geometric  mean 

life: 

49.3 

5.0  db  .  .  . 

changes  . 

•  • 

9.8 

Infei-red  slope 

of  fatigue  curve: 

a  =  3.9 
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SECTION  IV 


STUDIES  OF  PANEL  RESPONSE 


A.  PANEL  DESIGN  AND  MOUNTING 


A  structural  panel,  representative  of  flight  vehicle  construc¬ 
tion,  was  designed  and  a  number  of  specimens  were  made  by  the 
California  Division  of  the  Lockheed  Aircraft  Corporation.  A  drawing 
of  this  panel  Is  reproduced  In  Fig,  11.  Certain  somewhat  unusual 
aspects  of  the  design  were  dictated  by  the  available  siren  apparatus 
and  by  the  desire  to  eliminate.  In  so  far  as  possible,  complicating 
factors  such  as  effects  of  angle  of  Incidence  of  sound,  extremely 
nonlinear  response,  numerous  modes,  etc. 


The  panel  Is  made  by  riveting  together  two  0.020-ln.  sheets  of 
Ale lad  2024  aluminum  alloy.  The  flat  sheet  has  T-3  temper.  The  beads 
In  the  other  sheet  are  formed  when  the  sheet  Is  In  a  soft  condition 
and  the  piece  Is  then  heat  treated  to  T-42.  A  transverse  stiffening 
angle  of  l/l6-ln.  aluminum  (5/®  In.  high  by  3/^  In.  wide)  was  added 
at  midspan  to  the  design  shown  In  Fig,  11.  The  ends  of  the  assembled 
beaded  panel  are  riveted  to  strips  (0.125  x  2.00  In.)  of  aluminum. 

In  testing,  these  strips  are  bolted  securely  by  five  5/l6-ln.  bolts 
to  the  mounting  structure. 


Two  mountings  are  used,  one  for  pure-tone  excitation  eind  one  for 
random  noise  excitation.  For  pure-tone  tests,  each  of  the  aluminum 
strips  at  the  ends  of  the  panel  (Fig.  11 )  Is  bolted  to  a  sturdy  Iron 
angle  (2"x2"xl/4")  which  Is  In  turn  bolted  to  the  sides  of  the  1  ft 
by  1  ft  duct  down  which  sound  propagates  from  the  siren.*  The  Inside 
face  of  the  panel  Is  thus  exposed  directly  to  the  sound  pressure, 
which  propagates  In  a  direction  parallel  to  Its  surface.  The  outside 
face  (the  beaded  side)  Is  open  to  the  room  and  exposed  only  to  the 
relatively  low  level  of  ambient  noise. 

For  random  noise  experiments  the  aluminum  strips  at  the  ends  of 
the  panels  are  bolted  to  a  framework  of  Iron  angles  which  Is  secured 
to  the  open  top  of  a  plywood  box.  The  box  Is  made  of  one-inch  ply¬ 
wood  and  has  dimensions  l4"  x  21"  x  13".  The  purpose  of  the  box  Is 
to  protect  one  side  of  the  panel  from  sound  pressures.  In  laboratory 
measurements  a  single  box  Is  used.  In  the  field  measurements  behind 
a  Jet  engine,  three  boxes  are  secured  side  by  side  on  a  5  ft  by  11  ft 
cast  Iron  bed-plate.  Diagonal  baffles  of  l/2-ln.  plywood  are  arranged 
at  the  sides  of  the  boxes,  making  a  smooth  transition  to  the  bed-plate 
(see  Fig.  12).  It  Is  thought  that  the  baffles  cause  less  disturbance 
to  the  Jet  sound  field  than  the  boxes  alone;  no  experimental  compari¬ 
son  was  made . 


*The  siren  facility  used  In  these  tests  Is  described  In  detail  In 
Reference  31. 
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DESIGN  OF  THE  STRUCTURAL  PANEL 


FIG.  12  THE  PANELS  IN  THEIR  FIELD  MOUNTING 
FOR  JET  TESTS  (INSTRUMENTS  EX¬ 
POSED  FOR  THE  PHOTOGRAPH) 
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B.  PURE-TONE  RESPONSE 


The  response  to  pure-tone  excitation  of  the  panels,  mounted  in 
the  siren's  duct,  was  investigated  at  some  length  with  sound  from  a 
loudspeaker  in  the  duct  as  well  as  with  the  sound  generated  by  the 
siren.  The  motion  was  studied  by  eye,  with  stroboscopic  light,  and 
by  measurements  with  accelerometers  and  strain  gauges  in  many  places. 

It  is  found  from  low-level  measurements  that  there  are  two  prin¬ 
cipal  modes  of  response.  The  first  is  essentially  a  simple  beam 
motion,  the  panel  moving  in  and  out  of  the  duct  with  no  nodes  of  mo¬ 
tion.  The  second  is  a  twisting  motion  with  a  nodal  line  running 
lengthwise  down  the  center  of  the  panel.  Thus,  in  this  mode,  one  free 
edge  of  the  panel  moves  into  the  duct  while  the  other  free  edge  moves 
out . 


In  response  curves  taken  with  a  constant  low  level  of  sound 
pressure,  two  resonant  peaks  are  found.  One  peak,  typically  at  about 
110  cps,  corresponds  to  motion  principally  in  the  beam  mode.  The 
second  peak,  typically  at  about  120  cps,  corresponds  to  motion  prin¬ 
cipally  in  the  twisting  mode.  (A  third  peak  at  about  270  cps  corres¬ 
ponds  to  the  second  mode  of  the  beam  motion,  but  the  response  at  this 
frequency  is  very  much  smaller  than  the  response  in  the  other  modes.) 
Subsequent  Investigations  with  more  Intense  excitation  and  consider¬ 
ation  of  the  effects  of  the  abnormal  environment  in  the  siren's  duct 
(a  point  discussed  later)  indicated  that  the  lower  frequency,  beam¬ 
like  mode  was  the  one  of  principal  interest.  Therefore  the  response 
data  reported  below  all  pertain  to  that  mode. 

The  precise  frequencies  of  the  resonance  peak  and  the  relative 
proportions  of  beam  motion  and  twist  were  found  to  vary  considerably 
from  panel  to  panel.  For  example,  the  differences  in  the  strain  levels 
measured  on  the  top  and  the  bottom  beads  at  the  lowest  resonance  fre¬ 
quency  ranged  from  less  than  1  db  to  10  db  for  different  panels,  when 
measured  at  low  levels  (sound  pressure  levels  of  about  13O  db  re 
0.0002  mlcrobar) .  There  was  at  the  same  time  considerable  variation 
(about  +3  db)  in  the  ratio  of  the  sound  pressure  to  the  largest 
response  strain  from  panel  to  panel.* 

On  the  other  hand,  much  "cleaner"  response  was  found  in  all 
respects  at  higher  levels  of  excitation,  approaching  and  into  the 
fatigue  range.  Only  a  single  peak  of  response  was  found,  lying  at 
a  frequency  intermediate  between  the  frequencies  of  the  two  peaks 
found  at  low  level.  The  strains  measured  on  the  top  and  bottom 
beads  were  more  nearly  equal,  with  differences  in  level  generally 
less  than  3  db.  Moreover  the  ratio  of  sound  pressure  to  strain  was 
more  nearly  consistent  between  panels.  (In  one  notable  exception, 
panel  no.  5,  the  recorded  pressure-strain  ratio  was  5  db  smaller 
than  the  average  of  data  for  other  panels.) 


*The  largest  of  the  strains  in  any  of  the  beads  was  taken  as  the 
measure  of  panel  response,  since  this  study  is  concerned  with 
fatigue  lifetime  and  the  largest  strain  causes  the  earliest  failure. 
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Prolonged  investigation  of  the  responses  of  the  panels  at  high 
levels  of  excitation  would  have  destroyed  their  usefulness  for 
gathering  fatigue  data.  Moreover,  such  an  investigation  would  have 
been  a  diversion  from  the  purposes  of  this  study.  However,  many 
response  data  were  acquired  in  the  course  of  the  fatigue  tests. 

Response  curves  for  panels  numbered  17  through  20*  are  shown  in 
Fig.  13.  Each  curve  represents  the  average  of  a  large  number  of 
readings  (from  23  to  162)  taken  at  2  db  Intervals  of  sound  pressure 
level.  The  sound  pressure  was  measured  with  the  panel  blocked.  The 
nominal  strain  response  is  the  rms  of  the  signal  from  a  gauge  in  a 
standard  reference  position,  the  same  on  each  panel;  this  signal  was 
recorded  graphically  and  the  value  read  later.  Individual  readings 
made  on  the  same  panel  at  the  same  pressure  varied  over  a  range  of 
about  +  0.5  db  centered  on  the  average  value.  (In  a  few  Instances, 
particularly  at  lower  amplitudes,  variations  of  1.0  to  I.5  db 
occurred. ) 

It  appears  from  Fig.  13  that  the  strain  in  each  panel  is  pro¬ 
portional  to  pressure,  up  to  a  pressure  level  of  about  l46  db.  Beyond 
that  point  some  nonlinearity  appears.  However,  the  measurements  at 
the  higher  levels  are  not  as  reliable  as  those  at  lower  levels;  they 
were  made  during  very  short  Intervals  of  time.  Measurements  for  the 
highest  points  for  panels  19  and  20  were  made  in  2  seconds.  There  was 
not  time  for  careful  tuning  of  frequency  to  resonance.  In  view  of 
the  variations  of  the  frequency  of  peak  response  with  amplitude,  as 
noted  above,  there  exists  some  question  whether  the  measurements  were 
made  at  resonances.  Therefore  there  is  some  question  about  the  ac¬ 
curacy  of  the  indicated  amount  of  nonlinearity. 

The  response  data  of  Fig.  13  which  lie  in  a  range  of  sound  pres¬ 
sure  levels  from  l40  to  l46  db  is  summarized  in  Table  III.  We  propose 
to  use  the  average  difference  between  sound  pressure  level  and  strain 
level  as  a  "transfer  function"  in  comparing  various  predictions  of 
fatigue  life. 


TABLE  III 

PURE-TONE  RESPONSE  OF  PANELS  AT  RESONANCE 


Data  taken  from  Fig.  13 .  Pressure  levels  in  db 
re  0.0002  mlcrobar.  Strain  levels  in  db  re  10“°. 

Blocked  Sound  Pressure  Level 

Panel  Number  Less  Nominal  Strain  Level 

17 

93.3  db 

18 

93.6 

19 

93.0 

20 

94.1 

Average:  93.5  db 

*The  numbers  given  panels  are  for  identification  only;  they  do  not 
indicate  order  of  testing. 
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SOUND  PRESSURE  LEVEL,  PANEL  BLOCKED 
(DB  RE  0.0002  MICROBAR) 


152 


140  ' - 1 - 1 - 1 - 1 - L. 

47  49  51  53  55  57 

LEVEL  OF  NOMINAL  STRAIN  {DB  RE  10“®) 


FIG.  13  PURE -TONE  PANEL  RESPONSE  DATA 
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1.  Loss  Factor 


Measurements  and  calculations  were  made  to  determine  a  loss 
factor  for  the  panels  when  mounted  In  the  side  of  the  sound  duct 
which  leads  from  the  siren.  Casual  observation  indicated  that  the 
panel's  motion  was  well  damped  in  that  location. 

Determination  of  damping  from  the  bandwidth  of  the  resonance 
peak  was  made  difficult  by  the  close  proximity  of  the  two  resonances. 
The  peak  with  higher  resonance  frequency  (twisting  motion)  was  con¬ 
siderably  narrower  than  the  other  (besim  motion ) .  Prom  the  average  of 
measurements  made  on  a  number  of  panels,  the  loss  factor  was  estimated 
as 


pure-tone  loss  factor 

from  bandwidth  measurements:  =  O.O91  , 

for  the  beam  motion .  On  the  other  hand,  both  bandwidth  and  decay 
measurements  on  the  higher  mode  (twisting  motion)  indicated  that  its 
loss  factor  was  much  less,  about  O.O15. 

The  large  damping  for  beam  motion  and  the  much  smaller  damping 
for  twisting  motion  strongly  suggest  that  the  beam  motion  is  ab¬ 
normally  damped  by  virtue  of  its  environment,  i.e.,  that  acoustic 
radiation  deimping  dominates  internal  mechanical  losses.  An  Important 
corollary  is  that  the  response  In  the  beam  mode  will  Increase  marked- 
ly^  by  a  factor  of  about  6,  relative  to  that  in  the  twisting  mode 
when  the  panel  is  mounted  in  a  realistic  environment  without  the  ab¬ 
normal  radiation  loading. 

The  hypothesis  of  radiation  loading  can  be  checked  by  calcula¬ 
tions.  The  radiation  load  upon  the  panel  mounted  in  the  duct  con¬ 
sists  of  two  Impedances  in  parallel:  (l)  the  impedance  looking  down 
the  horn,  and  (2)  the  impedance  looking  back  up  the  duct  to  the  siren. 
The  first  can  be  estimated,  by  acoustic  horn  theory,  to  be  a  pure 
resistance  (a  specific  acoustic  resistance  equal  to  pc,  where  p  is 
air  density  and  c  is  the  speed  of  sound).  The  second  Impedance  is 
not  known  because  the  effective  Impedance  of  the  siren  source  has 
not  been  investigated. 

Suppose  that  the  impedance  looking  back  at  the  siren  is  reactive 
(no  power  flow).  Then  the  total  resistive  Impedance  presented  to  the 
panel  can  vary  from  zero,  when  this  shunting  reactance  Is  zero,  to 
the  value  computed  for  the  horn  loading  alone,  when  the  shunting  re¬ 
actance  is  large.  This  maximum  resistive  impedance  leads  to  maximum 
acoustic  radiation  damping.  Prom  the  dimensions  of  the  panel  and  an 
assumed  mode  shape  (fundamental  mode  of  a  pinned  beam)  we  have  com¬ 
puted  that  the  maximum  loss  factor  Is 

pure-tone  loss  factor 

due  to  acoustic  radiation:  ri^  ^  0.154 
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However,  suppose  that  the  acoustic  Impedance  looking  up  the  duct 
toward  the  siren  is  also  a  pure  resistance  equal  to  pc  (l.e.,  a  pro¬ 
gressive  wave  with  no  reflections).  Then  the  total  radiation  load 
will  equal  half  the  maximum  .load  of  the  preceding  paragraph.  The 
resultant  prediction  of  loss  factor  is  half  as  large: 

pure-tone  loss  factor 

with  (pc/2)  acoustic  loading:  ti.j.  =  0.077 

A  final  estimate  of  the  pure-tone  loss  factor  was  based  on  com¬ 
parisons  of  the  pure-tone  response  and  the  random  response  measured 
In  the  laboratory.  Details  are  given  in  the  next  subsection.  For 
comparison,  the  result  Is  given  here; 

pure-tone  loss  factor 

Inferred  from  random  response:  'H.j.  =  0.13 

The  four  estimates  for  pure-tone  loss  factor  have  a  range  which 
is  not  greatly  reassuring. 


C .  RANDOM  RESPONSE 

The  response  of  the  panels  to  random  noise  was  measured  in  an 
acoustic  environment  resembling  that  to  be  used  in  the  jet-noise  field 
tests.  These  measurements  serve  two  purposes.  First  they  determine 
the  loss  factor  which  is  to  be  expected  in  the  jet  environment. 
Secondly,  in  the  programmed-amplltude  equivalence  scheme  (subsection 
II-D),  they  determine  the  scale  factor  by  which  one  can  translate 
from  sound  pressures  used  in  the  fatigue  tests  In  the  siren's  duct 
to  the  equivalent  sound  pressures  of  jet  noise. 

Panels  were  mounted  In  one  of  the  boxes  described  above  in  sub¬ 
section  A  and  were  irradiated  with  steady  random  noise.  Measurements 
were  made  of  the  rms  value  of  the  nominal  strain  (l.e.,  the  strain  at 
a  standard  location  on  the  bead  having  the  largest  strain).  The  noise 
was  generated  by  exciting  a  loudspeaker  with  the  amplified  signal  from 
an  electrical  random  noise  generator,  after  filtering  to  a  one-octave 
bandwidth  (75  to  15O  cps). 

There  was  considerable  difficulty  in  the  way  of  a  good  experi¬ 
ment.  A  more  powerful  noise  source  would  be  required  in  order  to 
fulfill  the  ideal  requirements:  a  good  "slgnal-to-nolse"  ratio  in 
the  strain  signal;  acoustic  noise  signal  generated  without  peak 
clipping;  sound  pressure  uniform  over  the  panel's  surface;  and  no 
objects  In  the  close  vicinity  of  the  panel.  The  measurements  re¬ 
ported  here  were  made  in  a  compromise  situation  in  which  the  loud¬ 
speaker  was  six  Inches  from  the  panel  The  sound  pressure  level 
varied  by  several  decibels  over  the  face  of  the  panel.  The  results 
of  these  random  response  measurements  are  given  in  Table  IV. 
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TABLE  IV 


RANDOM  RESPONSE  OP  PANELS  (LABORATORY  MEASUREMENTS) 


Pressure  levels  In  db  re  0.0002  mlcrobar,  measured 
In  a  third-octave  band  centered  on  110  cps, 
panel  moving.  Strain  levels  In  db  re  10“° 


Panel  Third-Octave  Band  Pressure  Level 
Number  Less  Nominal  Strain  Level 


15 

87.5  db 

16 

88.5 

17 

88.3 

18 

90.7 

19 

89.9 

20 

Average : 

92.4 

89.5  db 

Corresponding  measurements  of  response  were  acquired  In  the  course 
of  the  field  tests  behind  a  Jet  engine.  The  results  of  analysis  of 
tape  recordings  are  shown  In  Table  V. 

TABLE  V 

RANDOM  RESPONSE  OP  PANELS  (PIELD  MEASUREMENTS) 


Pressure  levels  In  db  re  0.0002  mlcrobar,  measured 
In  a  third-octave  band  centered  on  110  cps,  With 
pemel  moving.  Strain  levels  In  db  re  10”° 

Panel  Third-Octave  Band  Pressure  Level 
Number  Less  Nominal  Strain  Level 


10 

11 

12 

13 


86.8  db 

87.1 

88.0 

89.0 

Average:  87.7  db 


The  difference  between  the  average  responses  found  In  these  two 
sets  of  measurements  may.  In  the  absence  of  definitive  Information, 
be  ascribed  to  any  one  of  a  number  of  causes:  differences  In  the 
sound  signals,  variations  to  be  expected  In  the  average  of  a  small 
number  of  measurements,  variability  between  Individual  panels,  dif¬ 
ferences  In  the  acoustic  radiation  loading  In  the  two  tests.  We 
favor  the  last  two  explanations.  In  any  case,  the  difference  between 
averages  and  the  variations  between  measurements  on  Individual  panels 
are  Indicative  of  the  precision  which  can  be  achieved  In  a  measurement 
which  Is  all-important  In  any  fatigue  life  prediction  scheme. 
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1.  Loss  Factor 


The  loss  factor  for  the  panel  in  the  laboratory  random  environ¬ 
ment  was  determined  by  measuring  the  rate  of  decay  from  steady-state 
resonant  vibrations  excited  by  a  pure  tone  sound.  Data  for  the  six 
panels  listed  In  Table  IV  result  In  an  average  value  of  loss  factor 

loss  factor  In 

laboratory  random  environment:  =  0.0442  +  0.0057  t 

where  the  number  appended  with  +  Is  the  maximum  range  of  variation . 

The  values  of  loss  factor  and  response  (Table  IV)  In  this  en¬ 
vironment  can  be  combined  with  the  pure-tone  response  measured  with 
the  panel  In  the  siren's  duct  (Table  III)  to  predict  a  loss  factor 
for  the  pure  tone  environment.  This  prediction  Is  based  on  Eq.  (I7b). 
The  spectral  density  of  sound  pressure  Is  computed  from  the  mean 
square  pressure  In  the  third-octave  band  and  the  bandwidth, 

The  predicted  loss  factor  Is 

loss  factor  In  pure  tone  environment  predicted 

from  response  data  and  loss  factor  In  laboratory 

random  environment:  =  0.13  i 

which  Is  In  fair  agreement  with  the  previous  estimates. 

2.  Statistics  of  Random  Response 

A  15-mlnute  recording  of  the  strain  response  to  random  noise 
excitation  In  the  laboratory  was  analyzed  on  a  Probability  Density 
Analyzer.  The  rms  strain  response  level  was  l6.6  db  re  10”°.  The 
results  of  the  probability  density  analysis  are  shown  In  Pig.  l4. 

They  are  seen  to  be  In  very  close  agreement  with  the  Gaussian  dis¬ 
tribution  curve  appiKjprlate  to  random  response  of  a  linear  one-mode 
structure. 

A  similar  analysis  was  made  of  the  strain  recorded  In  the  field 
during  the  Jet  noise  fatigue  tests.  In  this  case  the  rms  strain  level 
was  47.4  db.  The  probability  density  of  this  strain  signal  Is  com¬ 
pared  with  the  Ideal  Gaussian  distribution  In  Pig.  15.  There  Is  an 
Indication  of  some  response  nonlinearity,  the  density  being  very 
slightly  high  near  Iq  and  low  at  3a  and  beyond.  The  probability 
density  of  extrema  (strain  peaks),  being  the  negative  slope  of  the 
strain  signal  density,  appears  from  Pig.  15  to  be  very  slightly  higher 
than  the  Rayleigh  distribution  (the  negative  slope  of  the  Gaussian 
curve)  In  the,  region  above  about  1.5a.  These  results  are  In  quali¬ 
tative  agreement  with  the  type  of  "nonlinearity"  displayed  In  the 
pure-tone  response  curves  of  Pig.  13.  However,  the  differences  be¬ 
tween  these  results  and  the  Rayleigh  distribution  do  not  appear  very 
large  In  the  region  of  the  "peak  damage"  strain  which  Is  2.2a  In  the 
present  case  [Eq.  (2l),  taking  a=4  as  validated  below]. 
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PROBABILITY  DENSITY  OF  STRAIN  SIGNAL 


1.0 


0  icr  2cr  3cr  4cr  5C7  GO- 


AMPLITUDE  OF  STRAIN 


FIG. 14  PROBABILITY  DENSITY  OF  STRAIN  IN 

PANEL  EXPOSED  TO  RANDOM  NOISE  IN 
THE  LABORATORY 
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PROBABILITY  DENSITY  OF  STRAIN  SIGNAL 


FIG.  15  PROBABILITY  DENSITY  OF  STRAIN 
IN  PANEL  EXPOSED  TO  JET  NOISE 


ASD  TR  61  -  639 


52 


D.  COMMENTS  ON  NONLINEARITY 


In  these  response  studies,  with  siren  or  loudspeaker  and  with 
Jet  noise,  there  have  been  several  indications  of  nonlinearity. 

Among  these  are  the  statistics  of  the  response  to  Jet  noise  and  the 
curvature  of  the  pure-tone  response  curves  (Fig.  13)  acquired  in  the 
course  of  the  variable-amplitude  fatigue  tests.  The  latter  indica¬ 
tion  is,  as  noted  above,  somewhat  suspect  because  of  the  Impossibility 
of  tuning  accurately  to  resonance.  In  contradiction  to  it,  more  care¬ 
ful  but  less  numerous  measurements  made  in  some  of  the  constant- 
amplitude  fatigue  tests,  using  the  siren,  seemed  to  indicate  linearity 
at  high  levels  of  response.  Specifically  the  ratio  of  blocked  sound 
pressure  to  strain  was  constant  (+  0.3  db)  from  panel  to  panel,  with 
the  exception  of  one  rogue  panel  for  which  the  recorded  ratio  was  5  db 
lower  (panel  no.  5) • 

Other  indications  of  nonlinearity  are  found  at  lower  levels  of 
response,  as  recorded  above.  These  are  the  shift  in  resonance  fre¬ 
quency  and  the  merging  of  the  two  resonance  peaks  into  one,  as  level 
is  Increased.  The  shifts  in  resonance  frequency  were  much  greater 
than  can  be  explained  by  any  appeal  to  the  hard- spring  resonator 
type  of  response  appropriate  to  a  pinned  beam  (ref.  16). 

We  wish  to  postulate  a  qualitative  explanation  for  many  of  the 
observed  peculiarities,  in  support  of  which  there  is  not  sufficient 
evidence.  The  two  nearby  modes  of  response  found  at  low  level  should 
be  nonllnearly  coupled,  in  the  potential  energy,  by  virtue  of  their 
sharing  a  common  membrane  strain  and  tension.  With  Increased  exci¬ 
tation,  this  coupling  tends  to  modify  the  pure  tone  response  curve, 
making  the  two  resonance  peaks  approach  and  eventually  merge  because 
the  coupling  Increases.  It  is  quite  possible  that  changes  of  this 
sort  are  much  more  Important,  at  moderate  levels  with  nearby  reso¬ 
nances,  than  any  nonlinear  effect  on  the  strain  amplitude  at  the 
resonance  peaks.  Above  the  level  for  critical  coupling,  that  is,  for 
merging  of  the  peaks,  the  panel  would  behave  much  like  a  one-mode 
system,  until  at  very  high  levels,  nonlinearity  began  to  have  signi¬ 
ficant  effects  upon  the  amplitude  of  response.  The  validity  of  this 
postulated  explanation  should  be  investigated. 
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SECTION  V 


PROGRAMS  FOR  VARIABLE -AMP LI TUBE  FATIGUE  TESTING 


In  this  section  we  consider  the  Ideal  program  and  various  modi¬ 
fications  thereto  which  were  found  necessary  as  a  result  of  limitations 
In  the  existing  system  for  controlling  sound  pressure. 


A.  THE  "IDEAL"  PROGRAM 

As  explained  In  subsection  II-D,  It  Is  desired  that  the  sound 
pressure  amplitude  be  cycled  repetitively  through  a  "block"  In  such 
a  way  as  to  reproduce  the  actual  probability  density  of  the  response 
peaks  In  the  response  to  low-level  random  noise.  The  measurements 
reported  In  section  IV  showed  that  density  to  be  described  by  the 
Rayleigh  distribution.  We  desire  therefore  a  program  In  which  the 
pressure  envelope  Is  a  monotonlc  rising  or  falling  function  of  time 
and  has  a  density  given  by  the  Rayleigh  distribution.  Thus  the 
envelope  must  satisfy  the  differential  equation 

dt  =  +  (T/2)  X  exp[-x^/2]  dx  (22) 

where  x  =  p/a  Is  the  ratio  of  the  peak  amplitude  of  pressure  to  the 
rms  value,  and 


T  = 


dt 


The  solution  Is 

t  =  +  (t/2)  exp[-p^/20^]  ,  0  <  +  t  <  t  +  T/2  .  (23) 

In  the  Interval  of  time  -T/2  <  t  <  +  T/2,  the ^pressure  given  by  this 
equation  rises  from  zero  to  Infinity  (at  t=0)  and  then  falls  to  zero 
again  (Fig.  l6) .  This  defines  the  basic  block  of  the  program.  The 
duration  of  one  block  Is  T.  The  quantity  a  represents  the  long-time 
rms  value  of  the  pressure.  The  number  of  cycles  In  one  block  depends 
on  the  sound  frequency. 


The  slope  of  the  pressure  envelope  (Fig.  l6)  Is  Infinite  when¬ 
ever  the  pressure  Is  zero  or  Infinite.  Were  such  a  pressure  program 
used  to  drive  a  resonant  structure,  the  structural  response  would 
not  follow  the  pressure  In  a  quasl-steadystate  manner  when  the  rate 
of  change  of  pressure  amplitude  Is  very  large.  (We  assume  here  that 
the  sound  pressure  control  system  Is  Ideal  and  unlimited.)  Therefore 
the  program  must  be  truncated  at  a  high  level  (in  the  range  of  4a, 
say)  and  at  a  low  level  (at  some  small  fraction  of  la).  These  trun¬ 
cation  levels  must  be  chosen  so  that  the  omitted  peaks  would  not 
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FIG.  16  ENVELOPE  OF  PRESSURE  DURING  ONE  PROGRAM  BLOCK 
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have  made  a  significant  contribution  to  the  total  fatigue  damage. 

The  major  problem  appears  to  be  a  proper  choice  of  the  upper  truncation 
level . 

1 .  Distortion  from  Truncation  at  High  Levels 

An  analysis  has  been  made  of  the  distortion  of  the  response 
envelope  of  a  linear  simple  resonator  (one  mode)  to  an  excitation 
envelope  which  is  a  pyramidal  ramp  (Fig.  17).  The  distortion 
includes  (l)  a  time  delay,  which  is  not  of  importance,  and  (2)  a 
depression  of  the  maximum  response  amplitude  below  the  expected  max¬ 
imum  value.  The  amount  of  depression  depends  on  the  loss  factor  of 
the  resonant  system  and  the  rate  of  change  of  pressure  amplitude. 

(Our  analysis  assumed  a  constant  rate  of  change;  in  the  calculations 
we  have  used  the  maximum  rate  of  change,  as  calculated  from  Eq.  (22), 
which  is  a  function  of  the  truncation  level.)  For  small  values  of 
depression,  the  depression  in  decibels  is 

5  »  9-38  exp[x^/2]/x^'nN^  decibels,  (24) 

where  x  =  p/a  is  the  truncation  level,  q  is  the  loss  factor,  and  Nj 
is  the  number  of  cycles  in  one  block  of  the  program. 

It  can  be  seen  that,  for  the  depression  to  be  small,  either  the 
truncation  level  must  be  low,  or  the  number  of  cycles  in  a  block 
must  be  large.  Some  numerical  solutions  of  Eq.  (24)  are  given  in 
Table  VI,  for  a  loss  factor  ti  =  1/30,  a  typical  value  for  aircraft 
structure,  and  for  a  depression  6  of  0.5  decibels. 

TABLE  VI 

RELATION  BETWEEN  BLOCK  SIZE  AND  EXPECTED  MAXIMUM  RESPONSE 


Calculated  for  0.5  db  depression  of  maximum  with  loss  factor  t  =  I/30 


Expected  Maximum  Depressed  Maximum  Number  of  Cycles 

Response  Response  in  One  Block 

linear:  db  re  o:  (Depression=0 . 5db ) 


30 

9.5  db 

2.8a 

5.6x10 

3.5a 

10.9  db 

3.3a 

2.2x10' 

40 

12.0  db 

3.80 

l.lxio' 

50 

14.0  db 

4.7a 

6.1xl0' 

A  comparison  of  these  numbers  with  theoretical  distribution  curves 
for  damage  in  random  response  (Fig.  4)  indicates  that  it  may  be  neces¬ 
sary  to  omit  the  higher  response  peaks,  which  could  account  for  some 
20^  of  damage,  in  order  to  be  able  to  get  in  5  or  more  blocks  before 
the  structure  fails.  For  example,  with  2024  aluminum  alloy,  one  might 
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AMPLITUDE 


DELAY  TIME  :  T  =  Z  / 

oJq  =  RESONANCE  FREQUENCY 
r)  =  LOSS  FACTOR 

DEPRESSION  =  0.736  Xt 


FIG.  17  DISTORTION  OF  RESPONSE  ENVELOPE  NEAR 
THE  PEAK  OF  THE  PRESSURE  ENVELOPE 
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expect  a  total  life  of  0.5  to  1.0x10  cvcles  with  a  slope  of  the 
fatigue  curve  (at  the  peak  damage  point)  of  a  a  8.  Then  if  the 
pressure  envelope  is  truncated  at  4o  (with  a  maximum  response  of  3 *80), 
about  20^  of  deimage  will  be  omitted.  There  will  be  only  5  to  10 
blocks  in  the  total  fatigue  life.  In  order  to  get  higher  values  of 
response,  without  severe  distortion  of  the  distribution  of  peak  values, 
one  would  have  to  use  a  smaller  number  of  blocks .  But  a  large  number 
of  blocks  are  desirable  in  order  to  obtain  a  reasonable  representation 
of  ‘•he  random  response. 

The  considerations  of  the  preceding  paragraphs  illustrate  a 
fundamental  physical  limitation  in  the  attempt  to  use  programmed 
excitation  with  quasi-steadystate  response  of  a  highly  resonant 
structure.  The  only  alternatives  are  entirely  different  types  of 
testing:  either  (l)  excite  transient  response,  which  may  require  use 
of  a  Jet  engine  as  sound  source,  or  (2)  excite  non-resonant,  i.e. 
static,  response,  which  promises  formidable  experimental  difficulties 
with  a  large  structure.  Fortunately,  very  high  maximum  response 
levels  (e.g.,  50)  which  require  long  blocks  seem  only  to  be  required 
for  long-life  tests,  where  the  slope  of  the  fatigue  curve  is  small. 

2.  Acceleration  of  Testing  by  Truncation  at  Low  Level 

We  indicated  earlier  that  there  appears  to  be  no  difficulty  in 
achieving  controlled  resonant  response  to  values  small  compared  with 
the  rms  value.  However  there  is  the  possibility  of  saving  consider¬ 
able  time  and  achieving  an  accelerated  fatigue  test  by  truncating  at 
some  low  level.  The  only  problem  is  choosing  a  level  such  that  the 
effect  on  fatigue  life  will  be  small. 

If  fatigue  interactions  between  strains  of  different  amplitudes 
were  small,  the  choice  of  lower  truncation  level  is  readily  made  with 
reference  to  the  cumulative  damage  distribution  curves.  Fig.  4.  From 
those  curves,  we  see  that  only  a  few  percent  of  total  damage,  or  less, 
should  be  attributable  to  all  peaks  below  the  rms  value  for  slopes  of 
the  fatigue  curve  greater  than  a  =  4.  But  these  smaller  peaks  account 
for  39^  of  the  time  T.  [See  Eq.  (23).]  In  Table  VII  we  list  the 
truncation  levels  corresponding  to  10^  damage,  from  Fig.  4,  and  the 
corresponding  time  savings,  from  Eq.  (23),  for  a  range  of  values  of  a. 

TABLE  VII 

TIME  SAVED  BY  TRUNCATION  AT  10$^  DAMAGE  LEVEL 


Slope  of 

Fatigue  Curve, 
a 

Truncation 

Level 

(db  re  rms) 

Time 
Saved 
(Percent ) 

4 

+  3.4  db 

66.4^ 

8 

+  6.9 

91.5^ 

16 

+10.4 

99.6^ 
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Now,  the  choice  of  truncation  level  cannot  be  based  on  the  slope 
a  of  the  constant-amplitude  fatigue  curve  if  interactions  are 
importamt.  (if  interactions  are  not  Important,  there  is  no  need  for 
a  programmed  test.)  However  the  results  available  from  research  on 
fatigue  Interactions  (refs.  17,  25)  Indicate  that  predictions  based  on 
a  slope  of  4  with  vanishingly  small  endurance  limit  are  probably 
safe,  for  2024  aind  7075  aluminum  alloys.  Precisely  such  a  prediction 
is  entailed  in  the  first  line  of  Table  VII,  indicating  that  the  test 
may  be  run  in  one-third  of  the  lifetime  for  random  excitation,  by 
eliminating  all  peaks  below  1.5cr(+3.4db  re  a). 

It  is  quite  possible  that  even  more  accelerated  fatigue  tests 
would  be  reasonably  accurate  when  the  expected  lifetime  is  very  large. 
However  this  should  be  verified  experimentally,  by  performing  two 
tests  with  different  truncation  levels . 


B.  PRACTICAL  SIREN  CONTROL  PROBLEMS 

A  variety  of  means  was  investigated  in  the  attempt  to  achieve 
an  adequately  controlled  sound  pressure,  continuously  varying 
according  to  the  program  described  above.  These  tests  were  made  with 
the  24-ln.  diameter  siren  system  described  in  reference  31. 

Both  flexibility  aind  reliability  are  required  of  the  control 
system.  Flexibility  is  required  because  the  modulator  control  of 
the  siren  has  a  nonlinear  relation  between  input  electrical  signal 
and  sound  pressure  output,  and  that  relation  depends  upon  the  static 
air  pressure  in  the  supply  plenum.  Thus  the  control  signal  must 
be  tailored  to  fit  not  only  the  function  describing  the  desired  program 
but  also  this  nonlinear  relation;  the  adjustment  must  be  accomplished 
experimentally,  by  cut-and-try  methods. 

Reliability  is  required  of  the  system  because  the  rationale  of 
the  programmed  fatigue  test  is  based  on  the  repetitive  application 
of  the  same  basic  program  of  loading.  If  small  errors  were  detected 
after  running  one  block,  they  could  (with  suitable  flexibility  of 
control)  readily  be  corrected,  and  even  compensated  for,  before 
running  the  next  block.  However  correction  for  large  changes  would 
require  computation  based  on  a  fatigue  theory,  the  very  procedure 
which  programmed  testing  is  designed  to  avoid. 

Several  means  for  achieving  a  continuously  varying  sound 
pressure  were  Investigated.  An  electronic  cathode  ray  tube 
"Photoformer"  was  available  which  generates  an  electric  voltage  func¬ 
tion  proportional  to  the  shape  of  ajn  opaque  cut-out  placed  in  front  of 
the  cathode  ray  tube.  Awkwardness  in  adjusting  the  shape  of  the 
cut-out,  the  imprecision  in  such  adjustments,  and  long-time  d.c. 
drift  in  the  output  voltage  contraindicated  the  use  of  this  instru¬ 
ment.  An  Arbitrary  Function  Generator  (manufactured  by  Americam 
Measurement  +  Control  Company)  generates  an  electric  voltage 
proportional  to  the  height  of  a  curve  drawn  on  graph  paper  with 
conducting  ink.  Difficulties  in  achieving  reliable  signals. 
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changing  with  the  rapidity  required  near  the  peak  of  the  program, 
made  us  conclude  investigations  of  its  use. 

The  continuous  program  was  then  abandoned  in  favor  of  a  stepped 
approximation  thereto.  A  control  device  based  on  a  stepping  relay 
was  designed,  developed  and  proved  by  tests  with  the  siren.  This 
programmer  steps  up  and  down  a  total  of  l6  steps,  including  an  "off" 
position.  Both  the  dwell  time  and  the  level  of  the  output  signal  at 
each  step  are  adjustable.  Various  modes  of  operation  were  included; 
(l)  continuous  cyclingj;  (2)  stop  at  end  of  block;  (3)  stop  aind  turn 
off  at  end  of  step;  (4)  manual  control  of  time  of  stepping  (for 
setting  levels);  (5)  continuous  recycling  of  any  one  step  with  2 
sec  dwell  between  end  and  beginning  (for  setting  time  intervals). 
This  instrument  proved  to  have  the  requisite  flexibility  and  preci¬ 
sion.  With  it,  changes  in  sound  pressure  level  of  less  than  1  db 
per  step  can  readily  be  achieved. 

Unfortunately,  the  amplitude  stability  of  sound  pressure,  with 
the  siren  controls  fixed,  was  not  adequate.  This  represented  a 
deterioration  of  the  siren’s  performaince  from  its  original  state, 
apparently  due  to  difficulties  in  the  siren-modulator  control  system 
whose  source  was  not  readily  apparent.  There  were  short-time 
fluctuations  in  level  of  about  l/2  to  1  db  either  side  of  average 
and  long-time  drifts  of  about  1/2  to  1  db  in  an  interval  of  about  15 
minutes.  Occasional  surges  in  sound  pressure  have  occurred;  one 
surge  in  the  middle  of  a  fatigue  test  destroyed  the  specimen. 

Detailed  investigation  and  reworking  of  the  siren  control 
system  was  not  possible  in  the  course  of  this  project.  Therefore 
a  simplified  testing  program  was  devised,  with  manual  control  of  the 
sound  pressure  level  and  duration  for  each  step. 


C.  A  SIMPLIFIED  PROGRAM 

The  simplified  program  for  sound  pressure,  which  was  used  in  the 
programmed-amplitude  fatigue  tests,  was  designed  to  be  consonant  with 
the  capabilities  of  siren  systems  which  do  not  include  a  sound 
pressure  amplitude  modulating  device. 

The  basic  program  for  one  block,  designed  to  be  equivalent  to 
105  cycles  of  random  response  at  110  cps  (i.e.,  15  minutes),  consists 
of  8  minutes  of  excitation  in  nine  intervals  (Fig.  l8) .*  The  program 
was  applied  in  the  following  manner:  First  the  panel,  mounted  in  the 
side  of  the  duct  leading  from  the  siren,  was  blocked  from  vibrating 
by  moderate  pressure  from  a  diagonal,  sponge-coated  bar.  Then  the 


*  Please  note  the  two  ordinate  scales  in  the  figure.  The  "peak 

level"  is  the  logarithmic  measure  of  the  maximum  amplit\ide  of  the 
pressure  cycle.  The  "rms  level"  is  the  rms  value  of  the  pressure 
cycle  (3  db  smaller)  such  as  would  be  indicated  on  a  sound  level 
meter. 
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LEVEL  OF  PEAK  SOUND  PRESSURE  ,  DB  RE 


"blocked"  sound  pressure  level  was  set  to  the  appropriate  value 
specified  on  the  program.  The  blocking  bar  was  removed  for  the 
prescribed  time  interval,  and  quickly  replaced.  The  blocked  sound 
pressure  level  was  reset  to  the  appropriate  value  for  the  next  step, 
etc.  This  procedure  is  reasonably  straight-forward,  although  tedious. 
Suitable  experimental  details  for  its  easy  use,  depending  upon  the 
design  of  the  "siren  facility",  are  quickly  discovered  in  the 
application. 

The  time  scale  of  the  basic  program  can  be  scaled  to  fit  the 
expected  fatigue  life  and  the  desired  number  of  blocks;  that  is,  the 
duration  of  each  step  can  be  changed  by  the  same  factor.  There 
is,  of  course,  some  difficulty  in  maintaining  precision  in  the 
length  of  the  shortest  interval  for  a  short-life  test.  If  the 
expected  total  life  is  4x10^  cycles  and  4  blocks  are .desired,  the 
duration  of  the  highest  step  should  be  only  1.7  seconds. 

1 .  Design  Basis  for  Program 

The  program  shown  in  Fig.  l8  was  designed  in  this  manner. 

Previous  constant-amplitude  fatigue  tests  and  random  tests  in  a  Jet 
noise  environment  had  indicated  that  the  slope  of  the  fatigue  curve 
for  the  panels  (made  of  Alclad  2024  aluminum  alloy)  was  very  close 
to  a  =  4  in  the  range  of  interest.  (See  section  VI  below.)  The 
lower  truncation  level  was  then  chosen  to  be  a  peak  level  of  +1  db 
re  a;  this  corresponds  to  about  the  4  percentile  level  on  the  cumula¬ 
tive  damage  distribution  curve  for  a  =  4  (Fig.  4). 

After  several  tries,  balancing  the  impracticability  of  very 
short  intervals  and  of  level  changes  of  odd  size  against  the  desire 
for  moderately  small  changes  in  level  between  steps,  a  decision  was 
made  to  have  steps  differ  in  level  by  a  uniform  2  decibels.  The 
duration  of  each  step  was  immediately  determinable  [from  Eq.  (23)] 
since,  for  example,  the  first  step  must  last  for  the  time  spent  by 
the  ideal  program  in  the  interval  from  +1  to  +3  db  re  o.  However 
the  duration  of  the  highest  step  was  chosen  to  equal  the  total  time 
spent  by  the  ideal  program  above  9  db  re  a. 

The  level  of  each  step  except  the  highest  was  chosen  simply  as 
the  center  of  the  range  on  the  ideal  program  that  it  represents. 

Thus  the  first  step,  representing  the  Interval  from  +1  to  +3  db  re  o 
on  the  ideal  program,  is  run  at  a  peak  leveT  of  +2  db  re  o.  Any 
other  procedure  presumes  prior  knowledge  of  details  of  the  response 
and  fatigue  mechanisms  which  are  unavailable.  (Were  they  available, 
a  programmed  test  would  be  unnecessary.)  Moreover  any  adjustment 
for  a  presumed  daimage  distribution  represents  rather  picayunish 
fiddling,  since  the  adjustment  would  be  very  small  in  any  interval 
where  the  damage  density  is  large. 

Only  in  determining  the  level  of  the  highest  step  were  fatigue 
and  response  mechanisms  Invoked;  even  there  it  was  not  necessary. 

If  the  highest  step  had  been  designed  on  the  same  basis  as  ;he 
others,  to  represent  the  interval  +9  to  +11  db  re  o,  it  should  have 
had  a  level  of  +10.0  db  re  a  and  a  duration  of  15.5  sec.  Instead 
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the  duration  (17  sec)  was  chosen  to  represent  all  peaks  above  +9 
db  re  a,  and  the  level  chosen  to  represent  all  that  interval 
(assuming  linear  response.  Miner's  rule,  a  =  4) .  The  calculated 
level,  +10.04  db  re  o  was  rounded  off  for  use.  The  refinement  does 
not  make  a  great  difference  In  the  program,  as  a  practical  matter, 
and  In  theory  It  Is  probably  unjustifiable. 
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SECTION  VI 


PANEL  FATIGUE  TESTS 


Specimens  of  the  structural  panel  whose  design  and  response  were 
described  In  section  IV  have  been  subjected  to  sonic  excitation  until 
they  failed  In  fatigue.  These  tests  were  of  three  types:  Pure-tone, 
constant-amplitude  excitation  at  resonance  by  a  siren;  programmed 
variable-amplitude  excitation  at  resonance  by  a  siren;  aind  excitation 
by  Jet  noise. 

The  results  of  these  tests  are  given  In  this  section.  Detailed 
comparisons  aind  comments  on  the  results  are  given  In  section  VII,  to 
follow. 


A.  CONSTANT -AMPLITUDE  TESTS 

Seven  panels  were  tested  with  constant  sound  pressure  levels 
under  controlled  conditions.  (An  eighth  was  sacrificed  to  familiari¬ 
zation.)  Each  panel  cracked  cleainly  across  one  of  the  beads, 
approximately  at  midspan.  No  Interval  of  crack  growth  was  seen  In 
spite  of  continuous  observation  under  stroboscopic  light.  Each 
panel  except  number  2  cracked  on  am  exterior  bead,  that  Is,  near 
one  of  the  free  edges.  (The  crack  In  panel  2  passed  through  a  small 
dent  on  am  Interior  bead  although  prior  response  measurements  had 
Indicated  that  the  highest  nominal  strain  occurred  In  am  exterior 
bead.  This  behavior  may  Indicate  that  a  dent  In  a  bead  affords 
significant  strain  concentration.)  Except  for  panels  2  amd  6,  the 
crack  appeared  In  that  exterior  bead  which  had  been  shown  by  prior 
response  measurements  to  have  slightly  larger  strain. 

The  fatigue  lives  for  these  panels  are  plotted  In  Fig.  19  against 
both  the  peak  amplitude  of  strain  at  midspan  and  the  exciting  sound 
pressure  level  measured  with  the  panel  blocked. 

On  the  strain-life  plot  Is  superposed  a  fatigue  curve  for  Alclad 
2024  aluminum  alloy  taken  from  the  literature  (ref.  32).  (Values 
of  stress  were  converted  to  strain  by  a  value  of  Young's  modulus 
equal  to  1.06x10'  psl.)  It  was  noted  above  that  both  panels  2  and  6 
failed  at  unexpected  points.  Data  are  not  available  to  predict  the 
values  of  strain  for  these  points,  but  they  must  be  larger  than  the 
measured  values  used  In  plotting.  The  other  data  are  In  quite  good 
agreement  with  the  curve. 

The  pressure-life  plot,  by  comparison,  shows  more  scatter  of  the 
data.  Note  the  peculiar  results  for  panel  5^  whose  abnormally  lively 
response  was  noted  In  section  IV,  above.  Unfortunately,  no  reliable 
pressure  measurement  Is  available  for  panel  1. 

All  these  tests  were  run  with  unmodified  panels,  as  delivered 
by  the  manufacturer.  During  the  subsequent  field  tests  with  Jet 
noise.  It  was  found  expedient  to  drill  0.125-ln.  holes  In  the  center 
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of  each  bead  at  midspan.  Two  panels  with  such  holes  were  fatigued 
with  constant  sound  pressure  from  the  siren.  The  data  are  given  in 
Table  VIII,  along  with  the  value  of  strain  concentration  factor 
Inferred  from  comparison  of  these  results  with  the  smooth  curve  In 
Fig.  19. 


TABLE  VIII 

FATIGUE  DATA  FOR  STRUCTURAL  PANELS  WITH  l/8-IN.  HOLES 
IN  BEADS;  STRAIN  CONCENTRATION  FACTOR  FOR  HOLES 


Panel  10  Log  N  Nominal  Strain^evel* 

Number  (N  =  Life  in  Cycles)  (db  re  10-°) 


14  52.7  50.1 

15  51.9  50.1 

Average  52.3  50.1 

Conversion  to  level  of 

peak  strain  at  midspan*  +8. 2 

58.3 

Strain  level  for  failure  at  same  life, 
panel  without  hole  (from  curve  In 

Fig.  19)  66.5 

concentration  factor  8.2  db 

(in  linear  measure:  2.6) 


*Note:  The  "nominal  strain  level"  in  decibels  Is  the  logarithmic 

measure,  read  directly  from  the  face  of  a  calibrated  meter,  of 
the  rms  electrical  signal  from  a  bridge  contairfing  a  strain 
gauge  at  a  reference  point  on  one  of  the  beads.  This  reference 
point  was  not  at  midspan  because  gauges  at  midspan  frequently 
failed  before  the  panels  did.  However  the  conversion  factor 
was  repeated  checked  by  measurements  and  found  to  be  constant 
between  panels.  It  Includes  3-0  db  for  conversion  from  rms  to 
peak  level  for  a  sine  wave,  and  5*2  db  to  account  for  the  change 
In  position  of  the  gauge. 
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B.  PROGRAMMED  VARIABLE-AMPLITUDE  TESTS 

Four  structural  panels,  with  0.125-in.  holes  drilled  in  the 
centers  of  the  beads,  were  fatigued  according  to  the  simplified 
programmed  test  described  in  subsection  V-C,  above.  (One  failed  in 
the  process  of  such  a  test  due  to  a  malfunction  in  the  sound  pressure 
control  system.) 

The  use  of  the  program  given  in  Fig.  l8  starts  with  a  selection 
of  the  equivalent  long-time  rms  sound  pressure*,  that  is,  the  value 
of  a.  Then  the  duration  of  the  block  is  chosen,  based  on  the  expected 
lifetime  amd  the  Impractablllty  of  short  Intervals.  The  duration  of 
each  step  in  the  block  is  scaled  in  proportion.  The  pamel  is  mounted 
in  the  duct  of  the  siren  amd  exposed  to  the  prograimmed  levels  of 
blocked  sound  pressure  in  sequence,  step  by  step  as  previously 
described,  until  a  crack  is  observed. 

Panel  17  lasted  l8  blocks,  each  nominally  having  10^  cycles. 

Panel  l8  lasted  for  5  of  these  105-cycle  blocks.  Pamels  19  and  20 
were  tested  with  blocks  only  one-eighth  as. long  as  that  shown  in 
Fig.  l8,  that  is,  nominally  having  1.25x10^  cycles.  These  two  panels 
cracked  after  2-1/2  amd  3  of  these  short  blocks. 

1.  Correction  of  Raw  Data 


The  short  blocks  call  for  very  short  durations  at  the  higher 
steps  (2.1,  3-5j  7.1  sec).  Inevitably,  there  are  errors.  In  fact, 
the  graphic  records  of  response  show  that  the  average  times  spent  in 
the  higher  steps  were  too  long.  We  propose  to  correct  the  measured 
lifetimes  by  means  of  a  simple  linear,  non-interacting  fatigue 
theory.  The  theory  may  not  be  accurate,  but  it  should  be  accurate 
enough  since  the  corrections  are  small .  The  correction  formula  is 


N  =  N 


^actual^ 


a 


design  a 

^deslgn^ 


where  N 
N 

A 

X 

a 


design 


=  number  of  cycles  at  which  point  is  plotted 

=  number  of  blocks  times  nominal  number  of  cycles 
per  block 

=  total  time  (design  or  actual)  spent  on  step  x 
=  design  value  of  peak  amplitude  of  step 
=  slope  of  constant  amplitude  fatigue  curve. 


In  computation  we  took  a  =  4.  (See  Fig.  19) 


*  Because  of  the  stepped  approximation  and  the  truncation  of  small 
values,  the  value  of  a  is  not  the  actual  long-time  rms  sound 
pressure . 
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The  computed  correction  factors  on  life  are: 
panel  no.  19:  l.l6  (+0.6  in  10  log  N) 
panel  no.  20:  1.12  (+0.5  in  10  log  N) 

These  factors  were  applied  before  plotting  data. 

A  second  correction  factor  was  applied  to  all  data  to  correct 
for  small  variations  between  individual  panels  of  the  low-level 
resonance  frequency .  The  nominal  number  of  cycles  per  block  was 
computed  for  a  frequency  of  110  cps.  Actual  frequencies  ranged 
from  106  to  109  cps.  A  correction  equal  to  the  ratio  of  frequencies 
was  made  for  each  panel. 

Note  that  the  second  correction  would  not  be  necessary  if  the 
fatigue  life  were  measured  in  time  units  instead  of  cycles.  More¬ 
over  the  first  correction  would  then  have  been  different  (smaller). 
However  our  results  have  been  plotted  with  lifetime  in  cycles,  for 
easy  comparison  with  theoretical  predictions  and  other  data. 

2.  Results 


The  fatigue  lifetime  of  the  panels  subjected  to  programmed 
sotind  from  the  siren  are  plotted  in  Fig.  20  as  a  function  of  the 
equivalent  long-time  rms  sound  pressure  (measured  with  panel 
blocked).  The  results  have  a  surprising  consistency  and  show  little 
scatter,  in  distinction  to  the  constant-amplitude  results  plotted 
at  the  top  of  Fig.  19- 

The  curve  on  Fig.  20  labelled  "predicted  fatigue  curve"  requires 
some  explanation.  The  basic  fatigue  information  contained  in  the  curve 
is  the  constant  eimplltude  data  reported  elsewhere  in  the  literature 
(ref.  32).  From  this  we  constructed  a  strain-life  curve,  a  segment 
of  which  appears  in  the  plot  of  our  constant-amplitude  results. 

Fig.  19.  (The  reasons  for  using  strain  as  the  variable  in  place  of 
stress  were  explained  in  some  detail  in  subsection  II-A-2  and 
illustrated  in  section  III,  the  report  of  our  cantilever  beam  exper¬ 
iments.)  However  we  fotind  no  data  on  which  to  base  a  strain-life 
curve  for  strains  above  the  yield  point,  except  the  ultimate  strain 
which  causes  failure  in  one  cycle.  D* Amato  (ref.  7)  has  demonstrated 
the  continuity  of  the  strain-life  curve  (for  plain  2024  alimiinum)  and 
we  used  his  results  in  section  III  (e.g..  Fig.  8).  In  the  present 
case  (Alclad  2024)  interpolation  over  a  rather  wide  gap  was  necessary. 
That  part  of  the  resultant  curve  which  is  pertinent  to  these  results 
is  shown  in  Fig.  21. 

Also  shown  in  Fig.  21  is  a  prediction  from  the  constant- 
amplitude  data  of  a  random  fatigue  curve.  This  construction  was 
based  on  the  simple  theory  which  asstimes  a  Rayleigh  distribution  of 
peaks  and  the  nonexistence  of  fatigue  interactions. 

The  ordinate  for  the  random  fatigue  curve  is  the  level  of  the 
rms  value  of  strain  at  the  point  of  fatigue.  This  variable  must 
be  related  to  sound  pressure  level  before  the  prediction  can  be 
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BLOCKED  SOUND  PRESSURE 
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FATIGUE  CURVES  FOR  ALCLAD  2024  ALUMINUM 


compared  with  the  results  of  our  programmed  amplitude  fatigue 
tests.  Fig.  20.  A  constant  conversion  factor  was  used,  whose  value 
Is  based  on  the  measured  strain  concentration  factor  (Table  VIII ) 
and  the  measured  pure-tone  response  (Table  III) .  The  resulting 
prediction,  plotted  In  Fig.  20,  Is  in  quite  satisfactory  agreement 
with  the  results  of  the  programmed -amplitude  fatigue  tests. 


C.  JET  NOISE  TESTS 

Five  structural  panels,  with  0.125-in.  holes  drilled  In  the 
centers  of  the  beads,  were  fatigued  by  exposing  them  to  the  noise 
from  a  Jet  engine.  These  measurements  were  made  at  a  test  stand 
at  ASD  where  a  J75-P19  engine  was  mounted  horizontally  above  an 
extended,  flat  concrete  surface.  The  center  line  of  the  engine  was 
6  ft  above  the  concrete.  The  panel  mounting  was  described  in 

section  IV  and  is  illustrated  in  Fig.  12,  above. 

A  survey  of  the  sound  pressure  levels  was  made  in  the  region 
aft,  and  of  course  to  the  side,  of  the  engine.  Contours  of  constant 
sound  pressure  levels  in  a  third-octave  band  centered  on  110  cps 
(engine  at  military  power)  are  plotted  in  Fig.  22.  For  the  fatigue 

tests,  the  panels  were  placed  about  35  ft  aft  and  10  ft  to  the  side 

of  the  exit  nozzle.  Different  levels  of  excitation  were  achieved 
by  running  the  engine  at  different  speeds,  from  100  percent  to  about 
72  percent  of  full  speed  (rpm). 

In  a  test,  the  engine  speed  was  increased  until  the  chosen 
sound  pressure  level  was  reached.  If  any  change  in  the  pressure 
occurred,  the  engine  operator  was  instructed  to  adjust  speed 
accordingly.  Periodic  inspections  of  the  panels  were  attempted 
with  the  engine  running  at  speed,  but  only  inspections  with  the 
engine  idling  were  relied  upon  for  detection  of  cracks.  These 
were  made  at  intervals  varying  in  accordance  with  expected  lifetime. 

During  the  tests,  meter  readings  and  tape  recordings  were  made 
of  the  strain  and  acceleration  of  the  panels  and  of  the  sound  pressure 
level  3-in.  above  the  center  panel  (overall  level  and  level  in  a 
third-octave  band  centered  on  110  cps).  The  recordings  were  analyzed 
later  in  the  laboratory.  The  probability  density  analysis  of  the 
strain  was  shown  above  in  Fig.  15.  Some  response  data  (pressure/ 
strain  ratios)  were  given  in  Table  V. 

Narrow-band  (8^)  spectra  of  the  acceleration  at  the  panel's 
center  and  the  strain  on  a  bead  on  one  edge  are  shown  in  Figs.  23  and 
24.  The  acceleration  spectrum  resembles  very  much  the  expected 
result  for  a  one-mode  system.  However  the  height  of  the  resonance 
peak  cannot  be  used  to  determine  the  loss  factor  since  the  band¬ 
width  of  the  analyzer  (8^1  is  narrower  than  the  bandwidth  of  the 
structure  (estimated:  3^)*  The  strain  spectrum  shows  perturbations 
from  the  shape  for  one  mode.  The  bvunp  at  about  60  cps  probably  is 
due  to  hum  in  the  system. 
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AXIAL  DISTANCE  IN  FEET 
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CONTOURS  OF  CONSTANT  SOUND  PRESSURE  LEVEL  IN  A  THIRD- 
OCTAVE  BAND  CENTERED  ON  1 10  CPS  (J75-PI9  JET  ENGINE 
AT  MILITARY  POWER  ) 


LEVEL  OF  ACCELERATION  IN  8%  BAND,  DB  RE 


40  60  80  100  200  400 

FREQUENCY  IN  CPS 


FIG.  23  NARROW -BAND  SPECTRUM  OF  ACCELERATION 
RESPONSE  OF  PANEL  TO  JET  NOISE 
(145  DB  OVERALL) 
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FREQUENCY  IN  CPS 


FIG.  24  NARROW -BAND  SPECTRUM  OF  STRAIN 
RESPONSE  OF  PANEL  TO  JET  NOISE 
(145  DB  OVERALL) 
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The  bump  at  about  25O  cps  probably  reflects  response  In  the 
second  beam  mode,  whose  resonance  frequency  was  found  by  pure-tone 
measurements  to  be  about  270  cps  (section  IV). 

1.  Results 


Panel  no.  9  was  broken  in  an  unknown  time  less  than  two 
minutes.  Reliable  measurements  for  it  are  not  available;  the  data 
plotted  in  the  following  graphs  are  not  reliable,  representing 
attempts  to  use  some  fleeting  observations  of  meters.  Data  for 
the  other  four  panels  (lifetimes:  4.5,  27,  65,  and  65  minutes)  have 
been  taken  with  care  from  the  tape  recordings. 

Strain-lifetime  results  are  plotted  In  Fig.  25.  That  figure 
also  contains  a  "predicted  fatigue  curve"  which  Is  based  solely  on 
(1)  the  computed  random  fatigue  curve  In  Fig.  21,  and  (2)  the 
strain  concentration  factor  determined  by  pure-tone  siren  tests. 

Table  VIII.  Except  in  the  case  of  panel  no.  9,  the  agreement 
between  prediction  and  measurements  Is  quite  good. 

Sound  pressure-lifetime  results  are  plotted  In  Fig.  26.  The 
"predicted  fatigue  curve"  In  this  figure  Is  based  on  Fig.  21  and 
Table  VIII,  as  above,  and  upon  the  average  pressure-strain  ratio 
determined  in  the  field  (Table  V).  The  differences  between  data 
and  prediction  are  greater  here  than  In  Fig.  25.  (The  sound  pressure 
plotted  In  Fig.  26  were  measured  with  the  panel  vibrating,  as  were  the 
data  used  in  Table  V.  With  the  panel  blocked,  the  sound  pressure 
levels  In  a  third-octave  band  centered  on  the  resonance  averaged  1  db 
higher. ) 
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FIG.  25  STRAIN -LIFE  FATIGUE  DATA  WITH  JET  NOISE  EXCITATION 
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FIG.  26  SOUND  PRESSURE  -  LIFE  FATIGUE  DATA  WITH  JET  NOISE  EXCITATION 


SECTION  VII 


EVALUATIONS  AND  COMPARISONS 


When  the  evaluation  of  prediction  schemes  has  been  promised, 
there  comes  a  time  when  equivocation  over  the  detailed  significance 
of  data  must  be  set  behind  one  and  straightforward  predictions  com¬ 
pared  with  test  results. 

We  shall  compare  the  measured  fatigue  lives  of  panels  exposed  to 
various  levels  of  Jet  noise  with  predictions  obtained  by  two  schemes. 
The  first  prediction  scheme  is  that  described  in  subsection  II-D.  It 
is  based  on  measurements  of  fatigue  lifetime  under  excitation  by  a 
program  of  varying  sound  pressures  generated  by  a  siren.  The  equiva¬ 
lent  random  noise  level  for  a  specific  program  is  found  by  random-noise 
response  measurements  (ratio  of  rms  strain  to  rms  pressure.) 

The  second  prediction  scheme  is  based  on  computations  from  fatigue 
curves  taken  from  the  literature,  assuming  a  Rayleigh  distribution  of 
strain  peaks,  and  making  no  allowance  for  fatigue  interactions  between 
different  peaks.  Experimental  measurements  Include  (l)  constant- 
amplitude  fatigue  tests,  using  a  siren,  to  determine  a  strain  concen¬ 
tration  factor,  and  (2)  random-noise  response  measurements.  This  pre¬ 
diction  scheme  is  essentially  that  described  In  subsection  II-D,  except 
that  measurements  of  random-noise  response  are  used  in  place  of  compu¬ 
tations,  for  reasons  discussed  below. 

The  first  comparisons  are  shown  in  Fig.  27.  The  two  prediction 
schemes  give  very  nearly  the  same  results,  the  predicted  lifetimes 
being,  on  the  average,  long  by  a  factor  of  from  2  to  4  while  the  pre¬ 
dicted  sound  pressure  levels  for  a  given  life  are  high  by  about  2  db. 
The  random-noise  response  data  on  which  these  predictions  are  based 
are  the  results  of  laboratory  measurements  (Table  rv). 

When  the  predictions  are  based  on  random-noise  measurements  made 
in  the  same  environment  as  the  Jet-noise  fatigue  tests  (Table  V)  a 
different  picture  emerges  (Fig.  28).  Each  prediction  scheme  appears 
to  define  a  very  satisfactory  trend  curve  for  the  Jet-noise  fatigue 
data. 


The  dashed  curve  in  Fig.  28  Is  plotted  with  the  slope  determined 
by  fatigue  tests  of  randomly  excited  cantilever  beams  (Table  II).  It 
Is  evident  that  those  results  are  in  essential  agreement  with  results 
of  the  other  tests.  Therefore  these  cantilever  beam  tests  are  an 
acceptable  means  for  designing  accelerated  life  tests. 
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SOUND  PRESSURE  LEVEL  IN  THIRD-OCTAVE  BAND 
DB  RE  0.0002  MICROBAR 


10  LOG  N 


FIG.  27  JET  NOISE  FATIGUE  DATA  COMPARED  WITH 
PREDICTIONS  BASED  ON  LABORATORY 
MEASUREMENTS  OF  RANDOM  RESPONSE 
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FIG.  28  JET  NOISE  FATIGUE  DATA  COMPARED  WITH 
PREDICTIONS  BASED  ON  FIELD  MEASUREMENTS 
OF  RANDOM  RESPONSE 
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A.  GENERAL  COMMENTS 


We  present  here  a  number  of  comments  on  the  Implications  of 
these  results. 

The  extent  of  the  agreement  between  predictions  and  measure¬ 
ments  apparent  In  Fig.  28  must  be  considered  a  little  bit  fortuitous. 
The  predictions  are  based  on  a  series  of  measurements,  each  of  which 
Is  somewhat  Inaccurate  for  many  reasons.  The  number  of  samples  used 
In  both  laboratory  and  field  tests  Is  not  large  and  changes  of  0.5  to 
1  db,  or  perhaps  more,  might  be  expected  In  a  repeat  test.  Neverthe¬ 
less,  satisfactory  agreement  has  been  achieved. 

The  basic  problem  In  these  tests  was  In  achieving  a  satisfactory 
measurement  of  response.  For  example,  fatigue  curves  plotted  In  terms 
of  measured  strains  showed  excellent  agreement.  (See  Figs.  20  and  25, 
remembering  that  the  measurements  for  panel  9  are  suspect.)  To  some 
extent  this  problem  was  probably  due  to  variations  between  character¬ 
istics  of  Individual  panels.  It  was  very  much  accentuated  by  the  pres¬ 
ence  In  low  level  response  of  two  nearby  modes,  one  of  whose  Importance 
was  unrealistically  magnified  In  the  pure-tone  tests  by  the  abnormal 
acoustic  radiation  damping  of  the  other. 

A  particularly  Important  aspect  of  the  response  problem  Is  the 
determination  of  the  response  to  random  noise.  The  direct  experimental 
approach  In  the  laboratory  has  Its  difficulties  llabnormal  radiation 
loading,  non-uniform  pressure  distribution,  etc.)  discussed  In  the  sec¬ 
tion  on  response.  Measurement  In  a  suitable  field  environment  may  not 
always  be  feasible.  The  other  alternative  of  calculation  from  measured 
values  of  damping,  using  Eq.  (17),  appears  from  the  present  experience 
to  be  particularly  poor.  Loss  factors  could  not  be  determined  with  as 
much  accuracy  as  the  pressure-strain  ratio. 

It  Is  for  this  reason  that  a  computational  procedure  based  on 
measured  loss  factors  and  Eq.  (l?)  was  not  Included  In  the  predictions 
exhibited  above.  Four  determinations  were  made  of  the  loss  factor  for 
the  panel  In  the  pure-tone  environment  (Section  IV  above).  The  result¬ 
ing  predictions  for  fatigue  curves  would  have  ranged  from  1.2  db  lower 
than  the  predictions  In  Fig.  27  to  3»2  db  higher.  It  Is  felt  that  the 
procedure  adopted,  using  the  laboratory  measurement  of  random  pressure- 
strain  ratio,  was  the  most  accurate. 

Finally,  this  Investigation  Illustrates  the  advantages  In  making 
many  measurements  of  response  on  many  specimens.  Accumulated  errors 
of  2  to  4  db  or  more  could  easily  have  resulted  from  reliance  on 
particular  Individual  measurements. 
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B.  PROGRAMMED  AMPLITUDE  TEST  SCHEMES  -  PRO  AND  CON 


The  programmed  amplitude  test  schemes  proposed  and  tested  In  this 
study  have  various  advantages  and  disadvantages. 

The  simplified  program,  with  2  db  steps,  proved  feasible  and  gave 
good  results.  It  Is  a  good  way  to  achieve  an  average  pure  tone  response 
without  actually  making  measurements.  If  It  had  been  applied  blindly, 
without  so  much  study  of  response  at  a  low  level,  many  of  the  worries 
about  the  two  modes  of  response  would  have  been  avoided.  On  the  other 
hand.  It  requires  much  more  test  time  than  constant-amplitude  fatigue 
tests,  and  could  become  quite  a  chore  for  long  lifetimes. 

A  continuous  program,  such  as  Is  described  In  subsection  V-A,  was 
not  necessary  In  the  present  tests,  to  Judge  from  the  .results.  However 
with  reliable  automatic  control  of  sound  pressure  amplitude  and  of  fre¬ 
quency,  a  continuous  program  would  be  quite  a  bit  simpler.  Continuous 
and  automatic  control  would  probably  be  essential  for  determinations  of 
short  lifetimes. 


C .  FATIGUE  INTERACTIONS 

The  fatigue  measurements  made  on  the  panels,  made  of  Alclad  2024 
aluminum  alloy.  Indicate  that  fatigue  Interactions  were  not  a  signif¬ 
icant  factor.  Therefore  no  Judgement  can  be  made  of  the  effectiveness 
of  the  programmed  amplitude  tests  In  simulating  this  aspect  of  random 
fatigue. 

However,  the  results  of  the  cantilever  beam  tests  on  plain  2024 
aluminum  alloy  (section  III)  Indicate  a  difference  between  random 
measurements  and  predictions  vrhlch  may  be  ascrlbable  to  fatigue  Inter¬ 
actions  . 


D.  NONLINEARITY 

Nonlinearity  of  response  Is  not  believed  to  have  been  a  strong 
factor  In  the  present  tests,  except  In  respect  to  Its  mode  coupling 
effects.  There  Is  some  suggestion,  both  In  the  programmed-amplltude 
tests  and  In  the  Jet  tests,  that  nonlinearity  may  have  accounted  for 
a  flattening  of  about  1  db  at  the  high  levels  In  the  fatigue  curves 
of  Figs.  27  and  28.  However,  a  definitive  conclusion  Is  not  appro¬ 
priate  In  view  of  the  experimental  accuracy. 
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SECTION  VIII 


RESEARCH  ON  SIREN  FATIGUE  TESTING 


The  results  of  this  Investigation  Indicate  a  number  of  areas 
where  research  Is  needed  In  order  to  validate  and  Improve  the  reli¬ 
ability  of  fatigue  testing  with  siren  excitation.  ITiese  areas  are 
listed  here  with  brief  discussions  of  the  Information  that  Is 
desired. 

A.  EXPERIMENTAL  INVESTIGATION  OF  VARIABILITY 
FOR  TYPICAL  CONSTRUCTION 

It  Is  extremely  Important  that  queintltatlve  limits  be  set  on  the 
accuracy  required  of  a  testing  scheme.  Quite  probably,  these  limits 
will  be  set  by  variability  In  the  response  and  fatigue  characteristics 
for  nominally  Identical  samples  of  the  same  structure.  What  are  typi¬ 
cal  values  for  these  variations?  If  the  variations  seem  too  large, 
can  they  be  reduced.  In  practical  construction,  by  suitable  quality 
control  techniques? 

The  panels  used  In  this  Investigation  afford  one  example  of 
variation  In  the  low-level  response  characteristics.  Yet  this  varia¬ 
tion  seemed  to  disappear  at  larger  response  amplitudes  (section  rv). 

Thus  predictions  of  life  based  on  the  low-level  response  to  pure  tones 
would  have  been  much  more  Inaccurate  than  predictions  based  either  on  the 
high-level  response  or  on  the  rms  response  to  noise. 

How  much  variability  Is  there  In  the  relations  between  fatigue  life¬ 
time  and  nominal  response?  In  the  panels  used  In  this  Investigation,  the 
strain  concentrations  were  quite  uniform,  from  panel  to  panel,  and  the 
strain-lifetime  relations  show  little  scatter.  But  much  more  scatter  Is 
likely  If  the  fatigue  occurs  at  some  sharp  bend,  on  a  rivet  line,  etc. 


B.  EXPERIMENTAL  COMPARISONS  OF  CONSTANT  AMPLITUDE, 

PROGRAMMED  AMPLITUDE,  AND  RANDOM-RESONANT  FATIGUE  LIFETIMES 

The  tests  made  In  this  study  afford  some  comparisons  of  the  fatigue 
lifetimes  of  a  resonant  structure  excited  by  constant  amplitude,  pro¬ 
grammed  amplitude,  and  random  signals.  Many  more  such  comparisons  could 
be  desired.  Especially  valuable  would  be  very  carefully  controlled  com¬ 
parisons  made  with  quite  Idealized  specimens.  A  shaker-excited  resonant 
beam  Is  one  likely  design.  The  experiments  should  be  refined  to  Include 
controlled  force  of  excitation  and  controlled  nonllnearltles  of  response 
(both  In  stiffness  and  damping).  They  should  Include  careful  response 
studies  and  a  moderately  large  number  of  tests.  Including  some  with 
quite  long  lifetime.  Various  materials  should  be  tested. 
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More  experience  with  variable  aunplltude  fatigue  testing  of  struc¬ 
tures  would  be  desirable.  These  should  Include  tests  with  other  designs 
and  other  materials.  The  practicability  of  long-life  tests  should  be 
studied;  the  lower  truncation  level  of  the  progreim  can  be  much  higher  In 
such  tests  with  a  consequent  saving  of  time.  In  any  test  with  lifetime 
less  than  about  105  cycles.  It  appears  that  automatic  control  of  the 
progreim  and  of  the  frequency  should  be  Incorporated;  such  control  de¬ 
vices  should  be  developed  and  tested. 


C.  RANDOM  RESPONSE  OP  ONE-MODE  RESONATOR  WITH  NONLINEAR  DAMPING 

Available  results  from  this  and  other  Investigations  have  Indicated 
the  practical  Importance  of  nonlinear  deunplng  In  governing  the  response 
of  structures.  It  Is  a  disturbing  fact  that  there  Is  neither  experi¬ 
mental  nor  theoretical  Information  on  which  to  base  an  estimate  of  the 
effects  of  nonlinear  deimplng  on  fatigue  life.  Both  experimental  and 
theoretical  Investigations  should  be  started.  Information  on  the  prob¬ 
ability  density  of  response  peaJcs  Is  essential. 


D.  RESPONSE  OP  NONLINEARLY  COUPLED  PROXIMATE  MODES 

The  peculiar  behavior  of  the  panel  used  In  this  study  (section  IV) 
has  been  tentatively  explained  as  resulting  from  the  presence  of  two 
resonant  modes  with  closely  adjacent  resonance  frequencies  and  from 
their  nonlinear  coupling  through  the  potential  energy.  According  to 
this  hypothesis,  the  coupling  results  In  essentially  one-mode  response 
when  the  eimplltude  of  vibration  Is  large,  even  though  the  force-response 
relations  are  not  noticeably  nonlinear. 

Practical  structures  very  often  show  two  such  proximate  modes  In 
their  response  curves.  Theory  Indicates  that  they  should  be  nonlinear ly 
coupled.  Both  theoretical  and  experimental  studies  of  such  systems  are 
required  for  adequate  understanding  and  prediction  of  their  response. 

If  the  explanation  proposed  here  Is  right,  prediction  of  the  fatigue 
lifetime  of  such  a  system  would  be  greatly  simplified. 


E.  EXPERIMENTS  ON  PATIQUE  LIPE  OP  A 
STRUCTURE  WITH  TWO  DISTINCT  MODES 

Structures  with  two  distinct  resonant  modes,  well  separated  In 
resonance  frequencies,  are  also  common.  Much  discussion,  little  anal¬ 
ysis,  and  apparently  no  experimentation  have  been  applied  to  the  ques¬ 
tion  of  the  fatigue  life  of  such  a  structure.  Some  carefully  controlled 
experiments,  with  associated  analyses,  are  required  to  Investigate  this 
area.  Both  constant  eimplltude  and  random  excitation  should  be  used, 
with  different  levels  at  the  two  frequencies. 
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F.  FEW-MODE  MODELS  OF  PRACTICAL  STRUCTURES 


Reports  of  detailed,  careful  experimental  Investigations  of  the 
response  of  practical  structures  to  excitation  by  sound  are  seldom 
found.  Since  theories  based  on  the  presence  of  but  one  mode  are 
simple  and  available,  evidence  Is  often  found  to  confirm  the  domi¬ 
nance  of  a  single  mode  of  response.  The  evidence  Is  not  always 
sufficient:  A  single  peak  In  the  spectrum  of  response  to  random 
noise  could  Indicate  only  a  clumping  of  modes  In  a  limited  frequency 
region. 

Experiments  and  analyses  should  be  made  to  define  the  response 
characteristics  of  reasonably  simple  structures  In  considerable 
detail.  This  work  might  well  start  with  pure-tone  studies,  leading 
to  the  description  of  an  Idealized  (possibly  nonlinear)  model  with 
a  few  modes.  The  response  of  the  model  to  random  forces  could  then 
be  studied  by  analyses  and  analog  experiments. 


ASD  TR  61-639 


85 


REFERENCE  LIST 


1.  J.  W.  Miles,  J.  Aeronaut.  Scl.  753-762  (195^). 

2.  W.  D.  Dorland,  NASA  TN  D-611  (August,  I96I). 

3.  A.  K.  Head  and  P.  H.  Hooke,  Proc.,  Int'l.  Conf.  on  Fatigue, 
sponsored  by  Inst.  Mech.  Engrs.  (Gt.  Brit.)  and  Am.  Soc.  Mech. 
Engrs.,  November  1956,  pp.  301-304. 

4.  A.  L.  Karnesky  and  H.  C.  Schjelderup,  Shock  and  Vibration 
Bulletin  No.  25^  Part  II,  pp.  39-54,  Office  of  the  (U.S.) 
Secretary  of  Defense,  Washington,  D.  C.,  December  1957i 

A.  L.  Eshleman,  Jr.,  J.  VanDyke,  Jr.,  and  P.  Belcher,  ASME 
Paper  No.  59-AV-48,  presented  at  the  ASME  Aviation'  Conference, 

Los  Angeles,  March  1959;  M-  Belcher,  J.  D.  VanDyke,  Jr.,  and 
A.  L.  Eshleman,  Jr.,  Aero/Space  Engineering  I8  (No.  6),  24-30 
(June  1959);  H.  C.  Schjelderup,  pp.  19-25  ln~~i^coustlcal  Fatigue, 
ASTM  Special  Technical  Publication  No.  284  (American  Society  for 
Testing  Materials,  Philadelphia,  Pennsylvania,  i960). 

5.  R.  W.  Frallch,  NASA  Memo  4-12-59L  (June  1959). 

6.  K.  Gunn,  Aeronaut.  Q.,  ^  277-294,  (1955). 

7.  R.  D' Amato,  WADD  TR  60-175  (April  i960), 

8.  A.  M.  Preudenthal,  R.  A.  Heller,  and  P.  J.  O'Leary,  WADC  TN  55-273, 
Part  I  (June  1955). 

9.  R.  H.  Lyon,  "On  the  Vibration  Statistics  of  a  Randomly  Excited 
Hard  Spring  Oscillator.  Il",  J.  Acoust.  Soc.  Am.,  1395-1403 
(1961).  (Prepared  under  Contract  AF  33(6l6)-7l47. ) 

10.  I.  Dyer,  P.  W.  Smith,  Jr.,  C.  I.  Malme,  and  C.  M.  Gogos, 

ASD  TR  61-262  (March  I96I) . 

11.  Reference  10,  Section  III. 

12.  Reference  10,  Section  II. 

13.  S.  H.  Crandall,  Chap.  4  In  S.  H.  Crandall,  editor.  Random 
Vibration  (The  Technology  Press  of  the  Massachusetts  Institute 
of  Technology,  Cambridge,  Massachusetts,  1958). 

14.  Reference  10,  Fig.  44. 

15.  R.  W.  Hess,  L.  W.  Lassiter,  and  H.  H.  Hubbard,  NACA  RM  L55EI3C 
(July  1955);  L.  W.  Lassiter  and  R.  W.  Hess,  NACA  Report  I367 

(1958). 

16.  P.  W.  Smith,  Jr.,  J.  Acoust.  Soc.  Am.  33,  752-756  (1961);  Ibid, 
8i4-815  (1961).  (Prepared  under  Contract  AF  33(61d)-7147.") 


ASD  TR  61-639 


86 


REFERENCE  LIST  (CONT) 


17.  A.  M.  Freudenthal  and  R.  A.  Heller,  WADC  TR  58-69  Part  1 
(Jiane  I958). 

18.  R.  H.  I^on,  J.  Acoust.  Soc.  Am.  33,  815  (I96l). 

19.  Karl  Pearson,  Tables  of  the  Incomplete  F  -  Function,  (Cambridge 
University  Press,  Cambridge,  England,  19^6) . 

20.  P.  W.  Smith.  Jr.,  "Response  of  Nonlinear  Structures  to  Random 
Excitation,"  presented  at  SAE  National  Aeronautic  and  Space 
Engineering  and  Manufacturing  Meeting,  Los  Angeles,  12  Octo¬ 
ber  1961,  Paper  No.  430E. 

21.  P.  W.  Smith,  Jr.,  C.  I.  Malme,  and  C.  M.  Gogos,  J.  Acoust.  Soc. 
Am.  ^  1476-1482  (1961). 

22.  A.  M.  Freudenthal  and  R.  A.  Heller,  Chapter  8  in  Fatigue  in 
Aircraft  Structures,  edited  by  Freudenthal  (Academic  Press, 

New  York,  1956) • 

23.  A.  M.  Freudenthal  and  R.  A.  Heller,  WADC  TN  55-273,  Part  II, 

( October,  1956 ) . 

24.  A.  M.  Freudenthal  and  R.  A.  Heller,  WADC  TR  58-69,  Part  II, 

( J anuary ,  i960 ) . 

25.  R.  A.  Heller,  WADD  TR  60-752  (February,  1961). 

26.  H.  F.  Hardrath  and  E.  C.  Naumann,  pp.  125-138  In  Fatigue  of 
Mrcraft  Structures,  ASTM  Special  Technical  Publication  No.  274 
(American  Society  for  Testing  Materials,  Philadelphia,  i960). 

27.  R.  H.  I^on,  J.  Acoust.  Soc.  Am.  32,  716-719  (I960). 

28.  P.  J.  Roark,  Formulas  for  Stress  and  Strain,  (McGraw-Hill  Book 
Company,  New  York,  1954),  Table  XVII,  p.  344. 

29.  R.  E.  Peterson  In  WADC  TR  59-507  (August,  1959). 

30.  Reference  10,  Section  VIII. 

31.  G.  W.  Kamperman,  C.  H.  Allen,  J.  R.  Sharp,  and  L.  J.  Williams, 
WADD  TR  59-12  (December,  1959). 

32.  H.  N.  Cummings,  WADD  TR  60-42  (July,  I960),  Fig.  96. 


ASD  TR  61-639 


87 


o 


3?  PH 


m 


«  CL  ^  o 
'w  *-i  rf 


■a&^ 

5  •< 


O  gi 

pH  cd  m 

ii  n 

■^  **  [^  ** 


O  r- 


-  IJ 

U  _  “3  ^ 

§•33  a  Sg 

^  ^  dj 


& 


®  - 
Z  P^ 


5^4 


rH  hJ  H 


5:  >  ^ 


T 


If 

S£ 

■2&J 

A  ^ 

&  < 


rn 

s? 


^  a 

OQ  0 
0  O 
H  O 


^  ^  0)  tj 

0  e| 

-0  u  -  "S  ^ 

§■3318 

“  fl  9 

n 

s  S  ^ 

S  2  ft 


^  CQ 

Si" 

■<  n 

a  ■< 


►j  a 


-  “  F 

£ 

^5  -s 

■s  a^ 

—  Tj  3 

^  0  A 

ea  i*-i  H* 

s.a  I 


0  0 


I  N 

\  / 


no 

0 

0 

o 

lh 

n 

0 


«  0  tJ 

3  T?  S  5 


-3  ^ 

■H  “S 


/M  _ 

s  a 

a 
a 


^  t.  n 

a  -S 

^  s^*|b-S 

0  0  ^  0  fl  -Cl  3 

0  **-<^0  0 


II  § 


/  V 

\  / 


0 
0 
O 

tH  □ 

<D  0) 

O  -o 

a  I 

«  a 

«4H  0 


s  3 

*  0 

is 

(U 

■Q  5 
S;? 

s; 


K  a  I 

5  'S  0 

I  §  i 

^  o  I 
T3  ^ 

0  d 


0 

O 

*iH 

4-> 

i2 

0 

L< 

U 

O 

o 

0 

0 

o 


S  s 

1  s 

0  s 

^  s 


■0"  S 


i 


0  0 

ll 


2  B 


3  3 


I  0  g 

n  5  -S 


0 
Li  ^ 

O  0 


il 


o 

1 


